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Résumé : 


Dans cette thèse, nous avons mené une étude théorique des propriétés des semi-conducteurs 
magnétiques dilués de composition Mg,.,X,S (X =Cr, Cd, Zn, Co et Ni) dans la structure rock 
salt. Cette étude a été réalisée en utilisant la méthode des ondes planes augmentées linéarisées 
à potentiel complet avec les orbitales locales FP LAPW + lo. Cette approche repose sur la 


théorie de la fonctionnelle de la densité (DFT) mise en œuvre dans le code WIEN2K. 


L'approximation WC-GGA a été employée pour calculer les propriétés structurales et 
élastiques, tandis que l'approximation TB-mBJ (Becke-Johnson modifié par Tran-Blaha) a été 
adoptée pour évaluer les propriétés électroniques, magnétiques et optiques; quant aux 


propriétés thermoélectriques, elles ont été explorées en utilisant le code Boltztrap. 


L'étude structurale des composés a montré que les résultats obtenus pour le composé binaire 
concordent bien avec les résultats expérimentaux. Les caractéristiques électroniques et 
magnétiques des composés Mg,,X,S (x — 0.125, 0.25, 0.50, 0.75) indiquent qu'ils ont le 
caractère demi-métallique, avec une polarisation de spin de 100 % au niveau de Fermi et un 
moment magnétique total de 4 18. Tous les composés sont mécaniquement stables avec une 
forte anisotropie mécanique et leurs spectres optiques présentent une forte absorption dans la 


région du rayonnement ultraviolet. 


Les semi-conducteurs non magnétiques Mgo375X01259 (X — Cd, Zn) ont une large bande 


interdite et des propriétés optiques et thermoélectriques appréciables. 


Les spectres optiques des semi-conducteurs antiferromagnétiques de compositions 
Mgo875X0125S (X — Co, Ni) sont dominés par des pics de fortes intensités dans la région 


ultraviolette et présentent également des pics dans la région visible. 
Mots clés : 


FPLAPW+ Lo, WEIN2K, Semi-conducteurs magnétiques dilués (DMS), demi-métallicité, 


thermoélasticité, magnétisme, optique, TB-mBJ. 


Abstract: 


In this thesis, we conducted a theoretical study of the properties of diluted magnetic 
semiconductors with the composition Mg,.,X,S (X = Cr, Cd, Zn, Co, and Ni) in the rock salt 
structure. This study employed the full-potential linearized augmented plane wave method 
with local orbitals (FP LAPW + lo), grounded in density functional theory (DFT) and 
implemented through the WIEN2K code. 


The WC-GGA approximation was utilized to compute structural and elastic properties, while 
the TB-mBJ (Tran-Blaha modified Becke-Johnson) approximation was adopted for evaluating 
electronic, magnetic, and optical characteristics. Additionally, thermoelectric properties were 


explored using Boltztrap. 


The structural study of the compounds revealed that the results for the binary compound are 
in good agreement with experimental findings. The electronic and magnetic properties of the 
Me£.,X,S compounds (x = 0.125, 0.25, 0.50, 0.75) showed that they exhibit a half-metallic 
character with 100% spin polarization at the Fermi level and a total magnetic moment equal 
to 4 13. AÏl compounds are mechanically stable and exhibit strong mechanical anisotropy. The 
optical spectra of the different compounds show strong absorption in the ultraviolet radiation 


region. 


The non-magnetic semiconductors Mg, 375X91258 (X = Cd, Zn) have a wide bandgap and 


suitable optical and thermoelectric properties. 


The optical properties of antiferromagnetic semiconductors with compositions Mg, 535X0 1258 
(X = Co, Ni) are dominated by peaks of high intensity in the ultraviolet region with additional 
peaks in the visible region. 


Keywords: 


FPLAPW+ Lo, WEIN2K, Diluted Magnetic Semiconductors (DMS), ferromagnetic 
half-metallicity, thermoelasticity, magnetism, optics, TB-mBJ. 
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Introduction générale 


Introduction générale : 


Les semi-conducteurs magnétiques dilués (DMS) sont des matériaux intéressants car ils 
combinent à la fois des propriétés semi-conductrices et magnétiques. L'émergence des DMS a 
remarquablement impacté le monde de la spintronique. Dans les années 1980 à 1990, les 
chercheurs ont été orientés vers l’investigation de nouveaux matériaux par l’'introduction 
d’ions magnétiques dans des matrices semi-conductrices afin d’obtenir les propriétés 


électroniques et magnétiques dans un même matériau. 


Les DMS sont généralement composés d'un semi-conducteur hôte de type II-VI, IV-VI, II-V 
et d’un ou de plusieurs éléments parmi les métaux de transition ou de terres rares. Des travaux 
théoriques importants [1,2] ont débuté en 1945 avec l’avènement des semi-conducteurs dans 
l'industrie électronique. La découverte du transistor en 1947 par John Bardeen, Walter 
Brattain et William Shockley a révolutionné la technologie [3]. Les semi-conducteurs de type 
I-VI constituent une classe importante de matériaux car ils ont été à l’origine du 
développement de divers dispositifs électroniques et optoélectroniques. Lorsqu'ils sont dopés 
avec des métaux de transition comme le manganèse (Mn), le chrome (Cr) ou le cobalt (Co), 


ils donnent des propriétés intéressantes incluant le transport et le magnétisme. 


Les semi-conducteurs magnétiques dilués (DMS), obtenus par le dopage d’un 
semi-conducteur par des impuretés magnétiques tels que les métaux de transition [1,2] ou les 
terres rares [4], peuvent présenter une stabilité ferromagnétique ou antiferromagnétique, ils 
ont un double caractère (métallique et semi-conducteur) et possèdent une polarisation en spin 


complète au niveau de Fermi [2], ils sont aussi appelés matériaux demi-métalliques et leurs 


propriétés magnétiques et de transport en font d’eux de bons candidats pour les applications 


dans le domaine de la spintronique. 


La spintronique est une technologie fondée sur l’exploitation du spin des électrons et de leur 
charge électrique, les DMS trouvent de larges applications dans cette nouvelle branche de la 
physique. Les propriétés magnétiques découlant des DMS permettent de manipuler le spin, et 
permettent la réalisation de dispositifs plus performants et plus efficaces tels que les 
transistors spintroniques et les mémoires magnétiques [5]. La magnétorésistance géante 


(GMR) est la première application des DMS, après la magnétorésistance tunnel [6,7]. 


Les DMS sont envisagés pour être utilisés dans les capteurs magnétiques avancés [8] et 
possèdent des applications potentielles en optoélectronique [9], ils sont également étudiés 


pour leur potentiel dans les calculs quantiques [10]. 


A cause de l'intérêt fondamental des DMS, de nombreuses recherches théoriques et 
expérimentales ont été entreprises pour découvrir et développer de nouveaux matériaux 


appartenant à cette catégorie avec des caractéristiques plus performantes. 


Les méthodes théoriques ab-initio [11], basées sur la théorie de la Fonctionnelle de la Densité 
(Density Functional Theory - DFT) dans la physique du solide, permettent le calcul de la 
densité électronique et de l'énergie de l'état fondamental du système cristallin en résolvant 
l'équation de Schrôdinger [12]. Ces méthodes sont devenues des outils incontournables dans 
l’étude des systèmes complexes et ont connu un essor considérable car elles fournissent des 


résultats proches des résultats obtenus par l’expérience. 


Le travail dans le cadre de cette thèse porte sur l’étude de DMS réalisés par le dopage d’un 
semi-conducteur binaire appartenant à la famille II-VI par un élément magnétique, il est 


organisé comme suit : 
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Le premier chapitre comprend des généralités sur les semi-conducteurs du groupe II-VI et 
les semi-conducteurs magnétiques dilués (DMS), ainsi que leurs applications dans la 


spintronique. 


Le second chapitre est dédié à la théorie de la Fonctionnelle de la Densité (DFT) et aux 


approximations adoptées dans le calcul des propriétés des matériaux étudiés. 


Le troisième chapitre contient les résultats des propriétés structurales, électroniques, 
élastiques, optiques et magnétiques des matériaux Mg,,Cr,S (x= 0, 0.125, 0.25, 0.50, 0.75) 


dans la structure rock salt ferromagnétique. 


Dans le chapitre quatre sont présentés les résultats des calculs des propriétés 
optoélectroniques et mécaniques des semi-conducteurs Mg, 375X9125S (X= Cd, Zn) à large 


bande interdite dans la phase stable non magnétique. 


Le cinquième chapitre est dédié aux propriétés électroniques et optiques des semi-conducteurs 
Mgos875X0 1258 (X —Co, Ni) antiferromagnétiques. 

A la fin du manuscrit est présentée une conclusion générale résumant le contexte général dans 
lequel a été effectuée la recherche dans le cadre de cette thèse, les objectifs ciblés et les 


résultats réalisés ainsi que leur importance en termes d’applications. 
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CHAPITRE 1 
Les semi-conducteurs du groupe II-VT et les 


semi-conducteurs magnétiques diluëés DMS. 


Introduction : 


Au cours des dernières décennies, les semi-conducteurs du groupe II-VI du tableau périodique 
ont suscité un grand intérêt dans l'industrie moderne et ont été au centre de recherche des 
scientifiques [1,2,3,4] à cause de leurs importantes propriétés électroniques et optiques qui 
leur confèrent une large bande interdite et une large bandes d’absorption dans le domaine 
d’énergie couvrant la région du visible et de l’infrarouge. Le dopage de ces composés par un 
élément magnétique peut les transformer en matériaux dits semi-conducteurs magnétiques 
dilués (DMS) possédant un double caractère, métal et semi-conducteur. L'exploitation de la 
charge et du moment magnétique des porteurs de charges dans ces nouveaux matériaux 
conduit à un alignement parallèle ou opposé des moments de leurs porteurs spin- up ? et spin- 


down | leur polarisation de spins au niveau de Fermi est maximale. 


L'injection du spin dans les matériaux ferromagnétiques constitue la notion fondamentale de 
la spintronique qui est un domaine à larges applications entre autres, la magnétorésistance 


géante découverte en 1988, la magnétorésistance tunnel, spin-FET, etc. 


Aujourd'hui, la recherche est orientée vers l’investigation d’une nouvelle classe de DMS avec 


des propriétés améliorées pour de larges utilisations surtout en spintronique. 


L. 1. Généralité sur les semi-conducteurs du groupe II-VI: 


Les semi-conducteurs du groupe II-VI du tableau périodique suscitent un grand intérêt dans 
l'industrie moderne [1,2,3,4], ils sont formés par un métal du groupe 2 (ITA) ou 12 (IIB) (I = 


Me, Be, Cd, Hg ou Zn) du tableau périodique et un non-métal du groupe 16 (les chalcogènes, 
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VIB = O,S, Se ou Te). Les éléments du groupe II possèdent deux électrons de valence dans 
l'orbitale "s" la plus externe ainsi qu'une couche "d" qui peut être soit incomplète, soit 
complète, tandis que les éléments du groupe VI ont un total de six électrons de valence, 
répartis sur l'orbitale "s" ( deux électrons) et sur l'orbitale "p" (quatre électrons). Les éléments 
du groupe II donnent les deux électrons (cations) B°* (B=M£g, Be, Sr, Zn, Cd...) aux éléments 
du groupe VI (anion) A? (A=O, S, Se, Te) par l’hybridation sp3 des orbitales moléculaires. 
L'interaction entre ces orbitales conduit à une structure cristalline périodique, pouvant 
cristalliser dans la structure zinc blende (ZB), rock salt, ou wurtzite (WZ). Ces structures 
cristallines peuvent se former dans des conditions normales ou sous l’effet de contraintes telle 
que la pression et la température. Les composés Mes et MeTe présentent une stabilité dans la 
structure Rock Salt (RS)[5], par contre, les composés MgSe, CdTe, ZnSe, et ZnTe sont 


stables dans la structure cristalline de type zinc blend (ZB) [5,6]. 


Les semi-conducteurs du groupe II-VI offrent des propriétés électroniques et optiques très 
importantes et ont une large bande interdite allant jusqu’à 4 eV, ils sont utilisés en 
optoélectroniques dans les diodes électroluminescentes (LED), les diodes laser, les cellules 
solaires et les détecteurs optiques sensibles aux rayonnements des régions visible et 
infrarouge. Grâce à leurs excellentes propriétés, ces matériaux ouvrent la voie à des avancées 
technologiques prometteuses. Par exemple, les LED bleues/vertes sont à base de M£gsS ou de 
BaSe [2,3], tandis que les photodétecteurs infrarouges sont à base de MgTe ou de SrTe [1,7]. 
Toutes ces importantes applications découlant de l’exploration des semi-conducteurs du 
groupe II-VI expliquent les nombreuses recherches théoriques et expérimentales réalisées 


jusqu’à ce jour. 


Ce manuscrit est dédié au semi-conducteur binaire MgS appartenant à la famille des 


chalcogénures alcalino-terreux. Ce matériau a principalement été choisi pour ses propriétés 
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électroniques et optiques remarquables, qui seront présentées, étudiées et commentées le long 


du manuscrit. 


L.2. Les semi-conducteurs magnétiques dilués DMS: 


Les DMS sont connus avec le symbole A,-,M ,B, ils sont composés d'un semi-conducteur 


hôte AB appartenant aux groupes II-VI, IV-VL III-V, IV et II-IV-VI dopé avec des métaux de 
transition ou des terres rares qui sont injectés avec de faibles quantités [8]. L'introduction de 
l’élément M disposant d’une couche “d” incomplète (métaux de transition) ou d’une couche 
"4f" (terres rares) donne une variété de propriétés exceptionnelles au composé et ceci est 
principalement due aux interactions entre les moments magnétiques des porteurs localisés 
(élément M) et délocalisés (éléments de la matrice hôte) [9] (Figure. IL. 1.). Le premier DMS 
étudié est le composé Cd,,Mn,Te dans lequel des atomes de Mn” magnétiques occupent les 
sièges des atomes de Cd°* non magnétiques en faible quantité, ce composé cristallise dans la 


structure zinc-blende (ZB) comme le composé binaire CdTe [10] (Figure. I. 1.) 


Figure [.1. (a) Structure cristalline du semi-conducteur CdTe dans la phase zinc blende, (b) la 


substitution cationique de Cd‘? par Mn‘? dans la maille de CdTe. 
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Dans les DMS, les porteurs de charge (électrons ou trous) du semi-conducteur interagissent 
via l'interaction d'échange de Coulomb avec les spins localisés des ions magnétiques 
introduits [11]. Cette interaction d'échange induit un couplage entre les moments magnétiques 
localisés et les spins des porteurs itinérants. Elle conduit à un alignement parallèle ou opposé 
(antiparallèle) des spins des porteurs et des ions magnétiques selon que l'interaction d'échange 


d-d est respectivement ferromagnétique ou antiferromagnétique. 


Sous certaines conditions, notamment une concentration suffisante en porteurs de charge 
injectés optiquement ou électriquement, cette interaction d'échange peut donner lieu à des 
propriétés ferromagnétiques à température ambiante et à des températures supérieures à cette 
dernière, comme représenté sur la Figure I. 2. [12]. C'est l'impureté (dopant) qui est à l'origine 


du caractère "magnétique" des DMS. 


Le dopage de la matrice hôte (semi-conducteur binaire) par l'élément de transition M vise à 
produire des matériaux ferromagnétiques semi-métalliques (HM). Ces matériaux présentent 
un caractère métallique dans un canal de spin (up ou down) et un caractère semi-conducteur 
dans l'autre canal de spin. Les semi-métalliques sont caractérisés par un ordre 
ferromagnétique et par une polarisation de spin égale à 100% au niveau de Fermi [13]. Les 
différentes propriétés de ces matériaux en font d’eux de bons candidats pour des applications 


dans le domaine de la spintronique. 


L'ordre ferromagnétique dans les matériaux semi-métalliques est d'une importance capitale 
dans les applications technologiques, c'est pourquoi les recherches sont menées dans ce sens 
pour investiguer de nouvelles classes de ces matériaux et maitriser leurs diverses applications 


[14, 15, 16, 17, 18]. 


La découverte des semi-métaux a été faite par Groot dans les composés semi-Heusler 
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NiMnSb, PtMnSb [19] et dans les composés Heusler [20, 21, 22], elle a également été 
observée dans les matériaux à base de graphène [23], dans les matériaux pérovskite [13, 24], 
[25, 26, 27], dans les semi-conducteurs III-V [28,29] et II-VI [30,31] ayant subi un dopage 
avec des éléments magnétiques ou avec des terres rares, dans les oxydes alcalino-terreux [32] 


et dans les chalcogénures alcalino-terreux [18, 30, 33, 34, 35]. 


Le domaine des DMS a connu un intérêt croissant lors des dernières décennies en raison du 
potentiel applicatif de ces matériaux pour l'électronique de spin. Cependant, les mécanismes 
microscopiques à l'origine du ferromagnétisme dans les DMS font encore l'objet de 


recherches actives [11] 
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Figure I. 2. Schéma de trois types de semi-conducteurs : (a) les semi-conducteurs non magnétiques, 


(b) les semi-conducteurs magnétiques, (c) les DMS (ou semi-conducteurs de réseau hôte). 


L3. Classification des matériaux magnétiques: 


Les semi-conducteurs magnétiques peuvent être divisés en deux groupes distincts en fonction 
de leurs propriétés. D'une part, il y a les semi-conducteurs dans lesquels les éléments 
magnétiques s'organisent de manière ordonnée, comme les chalcogénures CdCr;,Se, [36 ,37], 


et certaines manganites; d'autre part, 1l y a les semi-conducteurs dans lesquels l’impureté 
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magnétique est aléatoirement distribuée. Les matériaux présentant le magnétisme peuvent être 
classés en différentes familles selon leurs caractéristiques d'aimantation du champ magnétique 


H et de la susceptibilité magnétique y,,, on distingue alors les cinq types suivants : 


L3.1. Le diamagnétisme: 


Le diamagnétisme est un phénomène magnétique fondamental observé dans la plupart des 
matériaux, bien que souvent très faible. Lorsqu'un matériau diamagnétique est exposé à un 
champ magnétique externe H, il développe une aimantation M opposée à la direction du 
champ magnétique, (Figure. IL. 3) créant ainsi une réponse magnétique négative. Cette 
propriété est le résultat d'une faible susceptibilité magnétique induite par le mouvement 
orbital des électrons dans les atomes constitutifs. En conséquence, le diamagnétisme induit un 
affaiblissement du champ magnétique dans le matériau, on note quelques matériaux 


diamagnétiques tel que le cuivre (Cu) et les non-métaux. 


Figure L.3. Arrangement diamagnétique des spins atomiques dans un cristal. Le moment magnétique 


de chaque atome (bleu) s’oppose au champ magnétique extérieur (violet). 


L3.2. Le paramagnétisme: 


Le paramagnétisme se manifeste lorsque les matériaux possèdent des moments magnétiques 
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atomiques non nuls. Sous l’effet d’une induction magnétique externe H, ces moments 
magnétiques s'alignent dans la direction du champ, provoquant une aimantation nette. 
Cependant, contrairement au ferromagnétisme, le paramagnétisme est une réponse temporaire 
qui disparaît lorsque le champ magnétique externe est retiré. Les matériaux paramagnétiques 
se caractérisent par une susceptibilité magnétique positive, ils sont fréquemment associés à 


des ions de transition non appariés tels que l'oxygène moléculaire et l'aluminium. 
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Figure I. 4. (a) Paramagnétisme des atomes libres, (b) arrangement paramagnétique des spins 


atomiques après application du champ magnétique. 


L.3.3. Le ferromagnétisme: 


Le ferromagnétisme est un phénomène au cours duquel une substance à la capacité d'acquérir 
une magnétisation permanente ou de répondre fortement à un champ magnétique appliqué AH. 
Le terme "ferromagnétisme" est lié au magnétisme, introduit par les éléments M contenant 
une “d” ou “f” partiellement remplies. Le fer (Fe, Fe‘*) est le premier élément où cette 
propriété a été identifiée, les nombres des spins-up et de spins-down sont différents[38]. 
Contrairement aux composés non magnétiques, les matériaux présentant le ferromagnétisme 
montrent un alignement parallèle spontané de dipôles permanents, résultant du mouvement 
des électrons autour des orbitales atomiques (Figure I. 5 (a)). C’est Pierre-Ernest Weiss qui a 


découvert les alignements structuraux de ces dipôles [39], bien que ces alignements entraînent 
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des magnétisations élevées même sans champ magnétique externe H, en dessous de 
température de Curie Tc. À mesure que l'on dépasse Te, les fluctuations thermiques 
conduisent à l'annulation de la magnétisation, car elles dépassent l'interaction entre les 
moments magnétiques résultant en un état paramagnétique comme c’est présenté dans la 


Figure I. 5. (b) 
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Figure L. 5. (a) Matériaux ferromagnétiques avec des moments magnétiques identiques et parallèles. 


(b)Transition de l’ordre ferromagnétique à l’ordre paramagnétique à la température de Curie Tc. 


Le matériau ferromagnétique peut se présenter sous deux formes selon la répartition 
électronique dans les bandes d'énergies tel qu'indiqué dans la Figure I. 6. Les matériaux 
ferromagnétiques et demi-métalliques ont respectivement une polarisation en spin inférieure 


ou égale à 50 % et inférieure ou égale à 100 %. 
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Figure IL. 6 Types des matériaux ferromagnétiques. 


L.3.4. L’antiferromagnétisme: 


L'antiferromagnétisme est une propriété qui se caractérise par un alignement opposé des 
moments magnétiques atomiques dans un matériau, ce qui donne une magnétisation globale 
nulle (voir la Figure 1.7). Cette interaction antiferromagnétique est principalement due au 
super-échange (J), un processus virtuel où un électron localisé fait "un saut" dans la bande de 
conduction, mais il est immédiatement remplacé par un électron de spin opposé provenant de 
la bande de valence. De plus, cette interaction dépend de la température et se manifeste en 


dessous de la température de Néel (T,) [40]. 
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Figure I. 7. Arrangement antiferromagnétique des spins atomiques dans un cristal. 
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L.3.5. Le ferrimagnétisme: 


Le ferrimagnétisme partage des similitudes avec l'antiferromagnétisme, sauf qu’il est 
représenté par deux sous-réseaux magnétiques ayant des moments magnétiques adjacents 
d’amplitudes différents. Cependant, au lieu de s'annuler mutuellement, ces moments 
magnétiques se combinent pour former un seul et puissant moment magnétique, parfaitement 
aligné sur le champ magnétique en environnant. Le ferrimagnétisme est fréquemment observé 
dans des composés présentant des structures cristallines complexes, tels que les oxydes 


métalliques. 


Lost 
+ 


Figure IL. 8. Arrangement ferrimagnétique des spins atomiques dans un cristal. 


1.4. Les métaux de transition: 


Les métaux de transition forment une catégorie essentielle d'éléments chimiques et occupent 
la partie centrale du tableau périodique, ils sont caractérisés par une couche « d » incomplète. 
Les métaux de transition appartiennent aux groupes IIB, IVB, VB, VIIB, VIIIB... jusqu'à IIB 
dans le tableau périodique (voir la Figure I. 9.), le remplissage progressif des couches ‘“d” 


donne une grande stabilité chimique. 


La fait que ces orbitales sont partiellement remplies, on obtient des propriétés magnétiques 


caractéristiques. Les éléments Cr, Mn, Co, … constituent une large famille de matériaux de 
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transition très étudiés, telles que les DMS. Les phénomènes les plus remarquables observés 
dans les DMS découlent de l'interaction sp-d entre les sous-systèmes magnétique et 


électronique. 


En raison de cette interaction, toute perturbation appliquée au sous-système magnétique, telle 
qu'un champ magnétique H où une variation de la température T, influe sur le sous-système 
électronique. Cela engendre une série de phénomènes magnéto-optiques et de 
magnéto-transport captivants, notamment une rotation de Faraday exceptionnelle, des 


polarons magnétiques liés, et des anomalies de magnétorésistance [41,42]. 
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Figure I. 9. Tableau périodique des éléments présentant les métaux de transition (en violet). 
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LS. Les interactions magnétiques: 


Les interactions magnétiques représentent le fondement des phénomènes magnétiques 
observés dans les matériaux. À l'échelle microscopique, elles émergent de l'interaction entre 
les moments magnétiques associés aux électrons au sein des atomes ; on distingue deux types 
d'interactions : les interactions d'échange de type sp-d et les interaction d'échanges de type d-d 
se produisant entre des ions magnétiques voisins, appelée aussi l'interaction d'échange 


d'Heisenberg et peut-être ferromagnétique ou antiferromagnétique. 


L5.1 Interaction d'échange entre porteurs localisés et délocalisés: 


Les porteurs délocalisés sont soit l'orbital ”s” du cation dans la bande de conduction, soit 
l'orbitale “p” de l'anion dans la bande de valence. Les porteurs localisés sont les orbitales “d” 
dans le cas des métaux de transition. L’interaction sp-d entre ces deux types de porteurs est 


définit par l’hamiltonien d’échange d'Heisenberg suivant (éq. I.1) : 


. Sp—d - . L . : “3% 
j 7": l'interaction d'échange entre les porteurs localisés et délocalisés. 


S; : est le spin des métaux magnétiques en position R;. 
s : est le spin des porteurs délocalisés en position r. 


Selon l'interaction entre les porteurs des bandes de valence et de conduction et tenant compte 
de l’approximation du cristal virtuel et de l’approximation du champ moyen, l’hamiltonien 


d’échange peut prendre les deux formes suivantes : 


HP =-xNB<S>Ss 
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H°714=-xNa<sS>s 


Où : 


x : la concentration de l'ion magnétique. 


Noÿ et Na : sont les constantes d'échange correspondant aux interactions d'échanges p-d et s-d 
respectivement, si Nof et Noa sont positifs l'interaction est ferromagnétique et si les signes 
sont opposés l'interaction est alors antiferromagnétique. En général, la valeur absolue de f est 


supérieure à celle de «. 


<$S > et s : sont les composantes parallèles du champ magnétique de spins des porteurs 


localisés et délocalisés respectivement. 


L5.2 L'interaction d'ions magnétiques: 


Les interactions magnétiques fondamentales identifiées dans les DMS sont : le super-échange, 


le double échange de Zener et l’interaction Ruderman-Kittel-Kasuya-Yoshida (RKKY). 


L5.2.1 Super-échange: 


Dans le cas de l'interaction du super-échange, la configuration électronique des ions non 
magnétiques est utilisée comme médiateur entre deux ions magnétiques voisins. Cela 
implique que les spins des porteurs localisés “3d” ne subissent pas d'interactions mutuelles 
directes, mais ces interactions sont assurées via un atome non magnétique à travers les 


porteurs “p” (O7, F, S? -). 


Ce mécanisme a été proposé par Kramer [43] et développé par Anderson [44]. Il décrit 
l'interaction de super-échange formalisée dans des matériaux ioniques tels que les fluorures 


[45]. Ces matériaux sont composés de deux ions magnétiques séparés par un ion non 
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magnétique. 


L'interaction  super-échange peut conduire à un couplage ferromagnétique ou 


antiferromagnétique selon le mode de distribution des états e, et t», des métaux de transition. 


L5.2.2 Double échange: 


L'interaction de double échange implique deux cations magnétiques, souvent des cations de 
métaux de transition, et un atome non magnétique. Les deux cations de valence participent à 
la formation de liaisons covalentes avec l’anion non magnétique. Cette interaction a été 
initialement proposée par Zener [46] pour expliquer le ferromagnétisme dans les composés 
LaMn;C0020;, et LaMnO:. La figure I. 10 (c) présentent un exemple du couplage double ou 
l'électron de l'orbitale” p ” de O7 partage l'orbital “d” de Mn“ créant une liaison covalente, et 
un électron de l'orbitale “d” de Mn*'partage l'orbital ” p” de l'oxygène O, selon la règle de 
Hund. Lorsque le spin d'un électron de l'ion Mn*'est aligné sur le spin de l'électron de l'ion 


Mn voisin, cela favorise le ferromagnétisme. 


Super Echange 


(a) 


+ 4 » 
The antiferromagnetic alignment Mn (d°) 0 (2p) Mn (d°) 


O (2p) Mn: (d?) 


Figure L.10. (a) Interaction Super-échange, (b) Recouvrement entre M-0-O-M, (c) Interaction double 
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échange : cas de deux ions Mn‘, Mn‘ [47,48] 


L5.2.3 RKKY (Ruderman-Kittel-Kasuya-Yosida): 


L'interaction de type RKKY (Ruderman-Kittel-Kasuya-Yosida) introduite pour expliquer le 
phénomène du ferromagnétisme des terres rares [49] décrit les interactions magnétiques 
indirectes entre les moments magnétiques localisés de l'orbital 4f de la terre rare et les 
électrons de conduction à travers un milieu non magnétique. Cette interaction, d'abord 
proposée par Ruderman, Kittel, Kasuya et Yosida dans les années 1950 [49], s’applique 


généralement à la magnétorésistance géante. 


La théorie du couplage RKKY induit des oscillations dans la force d'interaction magnétique 
en fonction de la distance entre les couches. Ces oscillations peuvent conduire à des 
alignements magnétiques alternés, favorisant soit des états ferromagnétiques ou 


antiferromagnétiques en fonction de la configuration géométrique. 


Concernant les semi-conducteurs appartenant au groupe III-V [50], le ferromagnétisme est 
expliqué par l'interaction RKKY. Lorsque la matrice semi-conductrice d'un semi-conducteur 
du groupe III-V est dopée par le manganèse, elle donne naissance à un semi-conducteur de 
type p. Cela induit un couplage d'échange entre les porteurs délocalisés (trous) présents dans 


la bande de valence et les ions Mn” localisés. 


L.6. Les DMS dans la spintronique: 


La présence du ferromagnétisme ou de l’antiferromagnétisme dans les semi-conducteurs 
magnétiques dilués leur confère des applications significatives dans plusieurs domaines, 
notamment les aimants temporaires à usage étendu en électronique, le développement de 


dispositifs avancés tels que les imprimantes tridimensionnelles et les ordinateurs haute 
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performance, les robots souples et les dispositifs électroniques souples et flexibles [51]. 


Ces matériaux ont également été efficacement utilisés dans le traitement du cancer par 
thérapie génique [52], dans les dispositifs d'intelligence artificielle [53] et dans le domaine de 


spintronique. 


L6.1 La spintronique: 


La spintronique, également appelée "transport électronique par le spin", exploite 
simultanément la charge et le spin de l'électron dans le but d'améliorer le traitement de 


l'information et la transmission des données. 


L'électronique est basée sur les propriétés essentielles des électrons [54] mais cet élément est 
aussi caractérisé par un spin dont la découverte a été faite en 1921 à travers l'expérience de 


Stern et Gerlach, permettant ainsi de quantification du moment magnétique d'un atome. 


La spintronique vise à tirer parti de la charge de l’électron et de son spin dans l’orientation des 
éléctrons et dans la production des courants. Nevill Mott [55] a prédi le rôle du spin dans le 
transport électronique. Cette notion a été reprise par Albert Fert et Peter Grünberg [56,57] qui 


ont découvert en 1988 la magnétorésistance géante (GMR) dans les systèmes Fe/Cr. 


Aujourd’hui, le spin est exploité intensivement pour donner de l’essor au développement de la 
spintronique dans beaucoup de directions, notament dans le domaine de la magnétorésistance 


tunnel, la spintronique moléculaire, etc. 


L6.2 Magnétorésistance géante (GMR): 
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La (GMR) décrit la variation significative de la résistance électrique observée dans une 
structure multicouche de matériaux magnétiques et non magnétiques sous l’effet d’un champ 
magnétique externe. La découverte de la GMR à la fin des années 1980 par Albert Fert et 
Peter Grünberg [56,57] a conduit à une percée révolutionnaire dans le domaine des capteurs 


magnétiques et des technologies de stockage de données. 


L'effet GMR se produit dans des matériaux se présentant sous forme de couches alternées de 
métaux ferromagnétiques séparées par des couches de métaux non magnétiques. Les couches 
ferromagnétiques ont des moments magnétiques qui peuvent être orientés parallèlement ou 
anti parallèlement l'un par rapport à l'autre. L'application d’un champ magnétique externe 
influence l'alignement des moments magnétiques dans les couches ferromagnétiques dans la 
configuration parallèle, ces moments s’alignent, par conséquent, parallèlement et les électrons 
peuvent circuler plus librement entre les couches ferromagnétiques, ce qui se traduit par une 
résistance électrique plus faible. À l'inverse, dans la configuration antiparallèle, où les 
moments sont antiparallèls, le flux d'électrons est entravé, ce qui se traduit par une résistance 


électrique plus élevée, comme illustrée dans la figure I.11. 


La sensibilité de la GMR à de faibles variations du champ magnétique en fait une technologie 
cruciale pour diverses applications, telles que les têtes de lecture magnétiques des disques 


durs et les capteurs de champ magnétique. 
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Figure IL. 11. Schéma illustrant le modèle électrique de la magnétorésistance géante (GMR). (a) 


Multicouches magnétiques Fe/Cr/Fe (b) Orientation parallèle et antiparallèle des couches.[56] 


L6.3 Magnétorésistance tunnel (MRT): 


La magnétorésistance tunnel (TMR) est constatée au sein des jonctions tunnel magnétiques 
(MTJ), des structures en couches composées de deux couches ferromagnétiques séparées par 
une fine barrière tunnel isolante. C'est dans ce type de matériau que l'effet Tunnel entre deux 


couches ferromagnétiques a été observé [58]. 


La principale caractéristique de la TMR est le changement significatif de la résistance 
électrique lorsque les moments magnétiques des couches ferromagnétiques sont orientés 
différemment. Lorsqu'une tension est appliquée aux bornes de couches ferromagnétiques, les 
électrons peuvent traverser la barrière isolante par effet mécanique quantique (l'effet Tunnel), 
de la couche libre à la couche de référence, et vice versa. La probabilité de l'effet tunnel 
dépend des états énergétiques des électrons dans les deux couches ferromagnétiques et dans la 
barrière isolante. Quand les moments magnétiques sont parallèles (configuration parallèle), la 


probabilité de l'effet tunnel est plus élevée, ce qui se traduit par une résistance électrique plus 
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faible. En revanche, quand les moments magnétiques sont antiparallèles (configuration 
antiparallèle), la probabilité de l'effet tunnel est plus faible, ce qui se traduit par une résistance 
électrique plus élevée. Ce changement de résistance est dû à l'orientation relative des 


moments magnétiques. 
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Figure I. 12. Schéma de la magnétorésistance à effet tunnel en configuration parallèle (haut) et 


antiparallèle (bas).[59] 


L.6.4 Spin-FET: 


Le spin-FET, également connu sous le nom de transistor à effet de champ spin (Spin Field 
Effect Transistor), a été proposé par Datta et Das en 1990 [60]. Il s'agit d'un dispositif de 
spintronique où le contrôle du flux de courant électrique est réalisé non seulement par la 
charge des électrons (comme dans les transistors à effet de champ traditionnels) mais aussi 
par leur spin intrinsèque. Le spin-FET est constitué d'un canal semi-conducteur, d'une source, 


d'un drain et d'une grille. Cependant, dans les spin-FET, l'orientation des spins électroniques 
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joue un rôle crucial dans la connaissance de la conductivité du canal. Le fonctionnement des 
spin-FET implique la manipulation de porteurs de charges polarisés en spin dans le canal 
semi-conducteur. Les porteurs polarisés en spin sont des électrons ayant une orientation 
préférentielle du spin. La tension de la grille est utilisée pour contrôler l'orientation du spin 
des porteurs dans le canal (la précession des spins d'une source ferromagnétique vers un drain 
ferromagnétique), influençant ainsi le courant électrique circulant entre la source et le drain, 


comme le montre la figure I. 13. 
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Figure I. 13 Schéma d'un transistor à effet de champ de spin.[60] 


L6.5 Spin-LED: 


Une diode électroluminescente de spin (spin-LED) exploite le spin des électrons pour 
contrôler et manipuler l'émission de lumière. Les diodes électroluminescentes (DEL) 
traditionnelles reposent sur la recombinaison des paires électron-trou pour émettre des 
photons. Dans les LED de spin, l'orientation des spins d'électrons est considérée comme un 


facteur supplémentaire influençant le processus d'émission. 
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Le principe de spin-LED est illustré par le schéma de la figure I. 14. Il est formé par une 
couche ferromagnétique (bleu pâle) pour injecter les électrons polarisés en spin dans la 
couche d'un semi-conducteur (orange) où se produit la recombinaison électron-trou dans la 
région active (jaune) pour produire de la lumière polarisée (rose, où la rangée indique la 


direction de la polarisation). 


Figure I. 14. Schéma de la SLED (jonction p-n) [61] 


L7. Sulfure de magnésium MsS: 


L7.1 Aperçu sur le MsS: 


Le sulfure de magnésium M£gsS est un chalcogénure alcalino-terreux [62] appartenant à la 
famille des semi-conducteurs du groupe II-VI [63], il existe en abondance dans la nature [64] 
et possède une large bande interdite, ce qui le rend intéressant d'un point de vue 
technologique. [65]. Dans des conditions normales de pression et de température, la structure 
la plus stable de Mg£sS est la structure rock salt (RS) [66,67] et sa bande interdite est indirecte 


[67]. Le MgsS peut également cristalliser dans la structure d'arséniure de nickel (NiAs) avec 
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une bande interdite indirecte, dans la structure de zinc blende (ZB) et dans la structure 
wurtzite (WZ) avec une bande interdite directe [68]. Le M£s est utilisé dans les dispositifs 
électroluminescents et optoélectroniques, en microélectronique, en catalyse, en médecine et 
en spintronique [69,62]. Le MgsS présente de faibles constantes diélectriques et des bandes 
interdites appropriées situées dans la région ultraviolette (UV) du spectre lumineux, ces 
caractéristiques sont recherchées dans les dispositifs optoélectroniques et photovoltaïques 


[70]. Drief et all. [S] ont étudié la structure la plus stable des composés MgS, MgSe et MeTe. 


Lorsque le M£gS est dopé avec des impuretés magnétiques, ses propriétés électroniques, 
élastiques, optiques et magnétiques sont modifiées [31,18,71]. Le dopage du MgsS par Cu‘? et 
Ce a démontré une bonne transition optique [72,73]. Des travaux sur le MgS dopé avec des 
métaux de transition (Cd, Zn, Cr, Ti, V, Ce.) [74,75,76,77,18,78] ont été réalisés pour 
l'investigation de nouveaux matériaux DMS [71,31], d’autres études théoriques portées sur les 
propriétés structurales et optiques des nanoparticules de M£S dopées avec Mn'° et Ag' [79,80] 


et l’étude du composé quaternaire Mg,Zn:,S,Se:., [81] ont également été réalisées. 


L7.2 Structure cristalline du MsS: 


Le sulfure de magnésium peut cristalliser dans cinq phases cristallines selon les conditions de 
pression et de température : Rock-Salt (B.), Chlorure de césium (B;) [82,83], Zinc blende 
(ZB) (B;), Wurtzite (WZ) (B;), et l'arséniure de nickel (NiAs)(B:1) [68]. Les études théoriques 
[5] et expérimentales [66] menées sur les composés MeTe, MeSe, MgS à la température 
ambiante et sous pression atmosphérique indiquent leur stabilité dans la phase rock salt FCC 
de groupe d’espace Fm3m N°225 avec une valeur du paramètre de maille égale à 5.198 À 
pour le MgS [84]. Dans cette structure, les cations (Mg?) occupent les positions (0,0,0) et 
l'anion (S?) occupe les positions octaédriques (2, 72,72.) dans la maille primitive (Figure. I. 


15 (c)); la maille conventionnelle possède 8 atomes : 4 atomes de Mg et 4 atomes deS, 
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comme le montre la Figure. I. 15 (a). 


(a) La maille conventionnelle 


Figure I. 15. Structures cristallines de MgsS dans la phase rock salt en format bille (a) (c), et en format 


polyédrique (b), obtenues par le logiciel VESTA. 


L7. 3 Structure de bandes du MsS: 


La configuration électronique des atomes Mg et S sont les suivantes : 

Mg : 1s°25°2p°3s7 

S : 1s22s22p°3s?p* 

La configuration des bandes électroniques de M£sS illustre la distribution des niveaux 
d'énergie au sein d'une bande. La bande de valence de ce composé est principalement 
constituée des orbitales 3s et 3p du soufre (S), ainsi que de l'orbital 3s du magnésium (M2). 
En ce qui concerne la bande de conduction, elle englobe les états électroniques d'énergie plus 


élevée de type 3p-$ et 3p-Mg. 
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L'écart énergétique séparant le haut de la bande de valence et le bas de la bande de conduction 
est appelé bande interdite. Le MgsS est un semiconducteur à gap indirect, car le minimum de 
sa bande de conduction (point 1°) et le maximum de sa bande de valence (point X) ne sont pas 
situés au même point (Figure I. 16). La valeur du gap énergétique du M£gs dans la phase RS a 
été calculée par Drief et all. [S] et par Zeyneb et all. [18], les valeurs respectives rencontrées 


par ces auteurs sont égales à 2.208 eV et 4.010 eV. 


Bande de Conduction 


BC 
E.(k) 


e_ 


Energie 


Bande de Valence 


Figure I. 16. Structure de bandes d’énergie d’un semi-conducteur à gap indirect. 

La structure des bandes est souvent représentée dans la première zone de Brillouin qui 
représente la maille élémentaire du réseau réciproque du cristal et qui est l’unité répétitive 
dans l'espace réciproque et dans laquelle sont contenues toutes les informations sur le réseau 
cristallin, les propriétés électroniques des matériaux sont généralement analysées dans cette 
zone. 

La limite de la zone de Brillouin est définie par des plans où les vecteurs d'onde ont des 
magnitudes spécifiques. La première zone de Brillouin de la structure rock salt est 
représentée dans la Figure I. 17. Les points de haute symétrie sont L' (0, 0, 0), X (1/2, 0, 1/2), 


W (1/4,1/4, 1/2), K (3/8, 3/8, 3/4) et L (0, 0, 1/2). Les lignes de haute symétrie sont À, À, et Z. 
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Figure IL. 17. Première zone de Brillouin de M£sS dans la structure rock salt. 
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Chapitre II 


Méthodes de calculs ab-initio basée sur la théorie de la 
fonctionnelle de la densité (DFT). 
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Introduction: 


Ce chapitre traite en détails les diverses techniques fondamentales utilisées dans le domaine 
de la modélisation moléculaire, il contient des aperçus approfondis de certaines méthodes et 
approximations utilisées dans la modélisation des systèmes quantiques complexes dont le 
traitement par les méthodes classiques de calculs s’avère très lourd et parfois impossible. 

Au début de ce chapitre est présentée l'équation de Schrôdinger (IL.1), équation fondamentale 
en mécanique quantique et dont la résolution passe inévitablement par des approximations 
adéquates et spécifiques aux systèmes traités pour l’obtention de solutions physiques. Ensuite 
est présentée l'approximation de Born-Oppenheimer (IL2), qui repose sur le principe qui 
consiste à séparer les mouvements nucléaires et électroniques pour rendre possible le 
traitement des systèmes moléculaires complexes. 

La deuxième approximation est celle d'Hartree et d'Hartree-Fock (11.3), fournissant un cadre 
conceptuel décrivant le mouvement des électrons sous l’effet d’un champ électrique moyen. 
Par la suite, la théorie de la Fonctionnelle de la Densité (DFT) (11.4), accompagnée des 
théorèmes fondamentaux d'Hohenberg et Kohn (1.4.1) et des équations de Kohn et Sham 
(1.4.2), toutes ces approximations ont été développées pour venir en aide à l’approche 
d'Hartree-Fock. 

Dans la théorie de la fonctionnelle de la densité, le seul problème qui entravait la solution de 
l’équation de Schrodinger était le terme lié à l’énergie d’échange-corrélation qui ne possédait 
pas une forme analytique pour la calculer. Dans ce contexte, un nombre assez important 
d’applications a vu le jour pour permettre d’arriver à une valeur de cette énergie et qui est 
tellement nécessaire pour résoudre cette équation, on peut citer, entre autres, les 
approximations de la Densité Locale (LDA) (1.5.1), du Gradient Généralisé (GGA) (IL.5.2), 
et le Potentiel de Becke-Johnson modifié par Tran-Blaha (TB-mbj) (ILS.5), offrant des 


compromis entre la précision et le coût du calcul. 


La fonction d’onde nécessaire pour la résolution de l’équation de Schrodinger est définie à 
travers des méthodes telles que la méthode des ondes planes augmentées et linéarisées 
(FP-LAPW) (IL6) ainsi que les méthodes APW (IL7), APW-+lo (ILS8), LAPW (IL9) et 
LAPW+LO (IL.10), chacune avec ses avantages spécifiques et ses limites d'application. 

Enfin, ce chapitre est clôturé par une présentation du code de calcul wien2k (11.12), un outil 


de calcul fondamental pour executer les méthodes de calculs qui offrent la possibilité d’une 
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compréhension plus large et plus approfondie des phénomènes quantiques complexes 


susceptibles de se produire dans les systèmes étudiés. 


IL. 1. Résolution de l'équation du Schrôdinger: 


L'équation de Schrôdinger, développée en 1927 par Erwin Schrôdinger, constitue le principe 
fondamental de la mécanique quantique, elle permet de décrire avec précision l'évolution de la 
fonction d'onde d'un système physique constitué de n électrons et de N noyaux à l’état 


fondamental. Cette équation, indépendante du temps, a la forme suivante [1] : 
Hyr, R) = Eyr, R). (IL.1) 


Où H est l'opérateur Hamiltonien exprimant l'énergie totale du système, ur, R) est 
la fonction d'onde électronique dépendant de la position r de l’électron et de la position R 


du noyau et E est l'énergie du système quantique. 

La résolution exacte de l'équation (IL.1) n'est envisageable que pour un système simple 
(mono-électronique), comme celui de l'atome d'hydrogène. Pour les systèmes poly 
électroniques (n électrons + N noyaux), il est impossible de résoudre analytiquement 
l'équation (II.1). 

L'hamiltonien H est composé de l'énergie cinétique des électrons a de masse m, et des 
noyaux Le de masse M; ainsi que de l’énergie potentielle due aux interactions 


coulombiennes répulsives entre électrons V,_, et aux interactions entre noyaux V et 
= nn 


attractives entre noyaux et électrons Vus FE s’écrit comme suit : 


à 2 
A 2 Vr 2 VR eZ Z 
h h ï 1 j 2 1 En) 
259 > m 2 > M ATE > sS BTE > F Le BTE > R°-R” 
i i 0 1] r, TR; 0 iÆj 0 iÆ) I J 


= = = = = (IL2) 
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Les méthodes basées sur des approximations dérivées de l'approche ab initio permettent de 
résoudre, de manière extrêmement précise, l'équation de Schrôdinger (11.1) et obtenir des 
solutions physiques plausibles. Suite à cela, il est possible de remonter aux structures des 
bandes électroniques des matériaux et, par conséquent, on obtient des résultats proches des 
ceux réalisés par des expérimentations. 

Dans la suite, nous allons présenter succinctement les approximations adoptées dans la 
résolution de l’équation de Schrodinger en commençant par l’approximation de 


Born-Oppenheimer jusqu’à arriver à la théorie de la fonctionnelle de la densité DFT. 


IL. 2. Approximation de Born-Oppenheimer: 


Dans les matériaux solides, les noyaux possèdent une masse très grande comparée à celles des 
électrons ; par conséquent, leur déplacement est beaucoup plus lent que celui des électrons 
dans le système. L'utilisation de l'approximation de Born-Oppenheimer (1927), simplifie 
efficacement la résolution de l'équation de Schrôdinger, elle repose sur l’hypothèse qui 
suppose que la position des noyaux atomiques est fixe par rapport à celle des électrons car ces 
derniers se déplacent avec une vitesse considérablement plus élevée. Dans de telles 


conditions, seule l’énergie cinétique des électrons T, est considérée, celle des noyaux Test 
par conséquent négligée et le potentiel répulsif entre les noyaux V_, est considéré comme 


constant [2]. 


Partant du fait que les mouvements des noyaux et des électrons dans une molécule sont 


séparés, la fonction d'onde |W)peut être exprimée sous la forme [3]. 


M) = (RluR()  (IL3) 


Où [Y (Rjket ly R(r) représentent respectivement les fonctions d’ondes nucléaire et 


électronique. 


L’Hamiltonien de l’équation (11.1), compte tenu de cette approximation, devient ce qu’on 


appelle l’Hamiltonien électronique ayant comme expression : 


“ h2 5 Tr, 1 5 eZ, 1 5 _ 1 5 22, a 4) 
Ron ne one ee gens A one Al one 2 Jes ur - 
e 2 : m, ATE, i A —* BTE, is) Ir. + | 8ne, dj [Re =R | 
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L'équation de Shrodinger s’écrit alors comme suit : 


Hi, =E V4, (IL.S) 
REA 1 eZ, 1 e 
En pe + One D DE N'ERODSEINUROI (LE 
Lo e BE SR, PSS 


L'approximation proposée par Born-Oppenheimer est restée insuffisante pour résoudre 
l'équation de Schrôdinger en raison de l'existence des interactions électron-électron, qui 
entravent la résolution totale du problème. Pour cela, il faut faire appel à d'autres 


approximations, telles que celles de Hartree et d'Hartree-Fock [4,5]. 


IL. 3. Approximation d'Hartree et d'Hartree-Fock: 


L'approximation d'Hartree (1928) [4] est essentielle dans la résolution de l'équation de 
Schrôdinger poly électronique en adoptant l’approche du champ moyen. Le mouvement de 
chaque électron est régi par le potentiel moyen engendré par les autres électrons et le noyau 
atomique. Ce phénomène est dû au fait que dans un système constitué de n électrons et de N 
noyaux, ils se comportent de manière indépendante, évoluant dans un champ électrostatique 


moyen sans être influencés par les spins ou les corrélations. 


Ensuite, la méthode d'approximation de Hartree-Fock (1930) est arrivée, [5,6]. Cette méthode 
améliore l'approximation de Hartree en prenant également en compte la symétrie des 
fonctions d'onde, ce qui permet d'introduire l'interaction d'échange. L'échange est causé par le 
principe d'exclusion de Pauli qui découle du principe d'Heisenberg. Celui-ci établit que toute 
fonction d'onde d'un système donné doit être antisymétrique par rapport à l'échange de deux 
particules. En 1930, Fock [5] suggère que la fonction d'onde d'un système composé de n 
électrons peut être définie par combinaison linéaire des fonctions d'onde des électrons 
individuels, tout en prenant la forme d'un déterminant de Slater qui respecte la règle 


d'exclusion de Pauli. 


d'Or) = WE) a? 
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AUN]AUPL AU) 
où -est le facteur de normalisation. 


et W :les fonctions d'ondes mono-électroniques. 


Les équations couplées du modèle de Hartree-Fock à résoudre deviennent : 


Fyf)-en(r) as 


L'opérateur Fock F : englobe l'opérateur cinétique, l'interaction avec le potentiel effectif, ainsi 


que le terme d'échange. 


Ces approches d'interaction de configuration ont été développées à partir de cette fondation. 
Le but est d'obtenir une solution exacte de l'équation de Schrôdinger. Toutefois, l'utilisation 
des calculs sur des échelles réduites est entravée par l'augmentation rapide du nombre de 
configurations lorsque le nombre d'électrons impliqués est important. Pour pallier ce 
problème. DFT a été développée. Elle facilite la résolution de l'équation de Schrôdinger en se 
servant de la densité électronique comme variable essentielle pour décrire l'état d'un système, 
à la place des fonctions d'onde. La simplification du système présente beaucoup d’avantages, 
car on arrive à réduire le nombre de degrés de liberté de 3N à ceux d'une fonction scalaire 
simple dans un espace tridimensionnel. Cela transforme efficacement un système à N corps en 


un système à une seule particule. 
IL. 4. Théorie de la fonctionnelle de la densité (DET): 


La DFT simplifie le problème quantique multi corps en utilisant un support indispensable 


pour simuler les propriétés d’un matériau. La principale idée de la DFT consiste à relier 


l'énergie totale d'un système électronique à sa densité électronique p (r) en utilisant une 
fonction d'énergie qui dépend de cette densité. Le développement de ce principe repose sur les 
travaux à l'échelle atomique de d'Hohenburg - Kohn et de Kohn-Sham [7,8] qui sont discutés 


dans la section (11.4.1) et (I1.4.2) 
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IL. 4.1. Théorèmes d'Hohenbersg - Kohn: 


L'approche de Pierre Hohenberg et Walter Kohn en 1964 [7] repose sur deux théorèmes 


essentiels : 


Théorème 1. « Pour chaque distribution de densité électronique p(r ) dans tout le système 


k : : ‘ ee : + : 
de particules, il existe un potentiel extérieur unique eus ) qui est associé aux noyaux 


dans leur état fondamental. » [7] 


Ce premier théorème démontre que l'énergie de l'état fondamental d'un système est 


entièrement déterminée par la densité électronique correspondante p{r ). Cela signifie qu'il 


n'y a aucune autre manière de décrire le système (fonction d'onde différente : Ÿ Æ  W) 


Théorème 2. « Le minimum de la fonctionnelle de l’énergie E rs minE p(r ) est donné par 


la densité de l’état fondamental p{r ), pour un potentiel externe v..(r) donné. » [7] 


Donc la fonctionnelle de l'énergie de la densité électronique obéit au principe variational 


selon l'expression suivante : 
ne —\ ,3 
E [pl=F,,lel +SV,.(r Jpfr dr, (11.9) 


Où la fonctionnelle de Hohenberg-Khon a est définie comme étant la somme de 


l'énergie cinétique (T) et de l’énergie potentielle des interactions coulombiennes répulsives 


entre les électrons ( Le ), elle s’écrit comme suit : 


FPT IP] +V,., le], (IL-10) 


e 


Selon ces deux théorèmes, l'approximation des équations de Kohn-Sham permet de 


Corriger et de reformuler l'erreur introduite dans l'expression de la fonctionnelle F uklP]- 
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II. 4.2. Équations de Kohn et Sham: 


Après les fondements de Hohenberg-Kohn, en 1965, Kohn et Sham [8] ont introduit le 
concept d'un système auxiliaire et fictif d'électrons non interactifs qui possède la même 
densité p(r) que le système réel. Les termes inconnus sont regroupés dans une fonctionnelle 


d'échange et de corrélation E XP} qui caractérise explicitement le système. Aïnsi, l'état 


fondamental à plusieurs corps peut être décrit avec précision. La fonctionnelle énergétique de 


l’état fondamental de l’équation (11.9) peut aussi bien s’écrire : 


p 
El Tpl+E,dPt Set (dr + L nv, | 11) 


= {EH[p]} 


OÙ 10 [ p] : est l'énergie cinétique classique des particules sans interaction. 

Vis (r) : L'énergie potentielle externe, 

E ALP] : L'énergie coulombienne classique de Hartree. 

E,, LP] : La fonction d'énergie d'échange-corrélation qui est une nouvelle fonctionnalité qui 


englobe tout ce qui n'est pas connu de manière exacte. En d'autres termes, elle est définie 
comme la différence entre l'énergie approximative et l'énergie exacte du système. Il est donc 


possible d'exprimer l'énergie d'échange-corrélation comme suit: 


cl 
ë, Lt [r tel 761] +", le - €, le] (1.12) 
Grâce à cette réécriture de l'énergie basée sur les N fonctions d'ondes mono électroniques, le 
problème devient plus simple et se réduit à trois équations, dénommées équations de Kohn et 


Sham [8] : 


1° équation de K-S : Définition du potentiel effectif 


PIC EP NI EM le Ne “‘(R) 


Avec : V, . [p (r” )] représente le potentiel d'échange-corrélation, donné par : 
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sE,e(r )] 


ARR ere 


(11.14) 


2°" Equation de K-S$ : introduction de V_.. [p (r)] dans les N équations de Schrôdinger 


Eff 


de fonction d’onde W: 


|- LV + Ve ŒO + V_ Ll + V,.lp ex |u, ()=E,W, (r) 


CV +7, lp GDU, Ge 4 (7) (LS) 


3°" Equation de K-S$ : relation entre la densité électronique et les N fonctions d’onde W. 


N 
p ()= 24") 2 (1.16) 


La résolution des équations (11.13), (IL.15) et (11.16) de K-S devient possible si l'énergie E Le 
[p] et le potentiel Fe à [p (r” )] d'échange-corrélation sont connus. 


La méthode itérative est utilisée pour obtenir le potentiel V, ff [p (r) à partir d'une densité 
électronique initiale (P,). En résolvant l'équation de Kohn-Sham (11.16), on calcule une 
nouvelle densité (p,) de USE Si cette dernière diffère de (P; ), les deux densités ( (P; jet (p, 


) sont combinées pour obtenir une nouvelle densité (P, ). L'exemple donné dans la Figure. II 


1. Illustre le schéma du cycle SCF utilisé pour résoudre les équations KS. 
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Résolution de l'équation(- 112 V2+ 
VEff{p(r)] ) &i(r) =ei Oi(r) 


Calcule de la nouvelle 
densité p1 


Calcule de VEff[p(r)] 


Calcule p2= pi+1 
pi 
Densité électronique 
initiale pi Convergence 


Figure II. 1. Représentation schématique du cycle SCF pour résoudre les équations KS 


Commencer avec 


IL. 5. Les approximations du potentiel d'échange-corrélation: 
L’obtention de solutions aux équations de Kohn-Sham est dépendante du terme d'échange et 
de corrélation E [p], qui reste l'inconnue et dont la connaissance facilite la solution de ces 


équations. Malheureusement, cette limite majeure de la théorie de la fonctionnelle de la 
densité n'a pas été clairement indiquée dans l'équation IL.12. L’absence d’une expression 


analytique pour Va conduit les chercheurs à proposer des modèles pour calculer d’une 


manière approximative ce terme sans lequel des solutions ne seraient pas possibles. De 
nombreuses approximations ont vu le jour, entre autres on cite l'approximation de la densité 
locale (LDA) et l'approximation de la densité de gradient généralisé (GGA). Dans la suite, 


nous discuterons dans les détails ces approximations. 


IL. 5. 1. L’approximation de la densité locale (LDA): 


La LDA, proposée par Kohn et Sham [8], est essentielle pour trouver une expression 
analytique satisfaisante du potentiel d'échange-corrélation tr). Cette approche 


suppose que le potentiel d'échange-corrélation est une fonction de la densité électronique 
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locale. On assimile le système à un gaz homogène d’électrons de densité p (r” ) en un point 


r. En utilisant cette approximation, on peut postuler que la fonctionnelle 


d'échange-corrélation prend la forme suivante : 


Ee EP = fp(r ke [P(r)d'r, (1.17) 


; DD de Pins Re 2 | 
Noté que £,. est l'énergie d'échange-corrélation du gaz d’électrons homogène, construite par 


ÿz LDA : À ; ; De 
deux termes, le terme d'échange E, P (r_) qui peut être calculée analytiquement en utilisant 
la fonction d'énergie d'échange de Dirac [9] et le deuxième terme est le terme de corrélation 
LDA ; FE ; 
EP (r ) qui est calculé numériquement par des méthodes de Monte-Carlo, par Ceperley et 


Alder [10]. 

L'estimation de la densité locale est particulièrement précise pour les systèmes covalents et les 
métaux simples, où la densité électronique varie lentement. Cependant, il faut souligner 
quelques inconvénients notables, tels que la sous-estimation du paramètre de la maille, de 
l'énergie du gap et des moments orbitaux de spin [11]. Par conséquent, il est essentiel d'ajouter 
des termes supplémentaires dans l'expression d'échange et de corrélation pour pallier ces 


problèmes, ce qui sera discuté dans la section suivante. 


IL 5.2. L’approximation du gradient généralisé (GGA): 


La GGA (Generalized Gradient Approximation) consiste en écrire la fonctionnelle d'échange 


À GGA : = ; L2t 
et de corrélation ss [p (r )] de la même manière que la LDA, mais avec une densité 


d'énergie d'échange et de corrélation p(r ) qui dépend à la fois de la densité électronique et 


du module de son gradient (pr ) ) [12]. Dans ces conditions, son expression est la suivante 


GGA 3 
EX Le œ-{b(r)vp{r)ldr , GLS) 
En général, la fonctionnelle d'échange-corrélation des systèmes polarisés en spin a la forme 


suivante : 


52 


EC Let@ pt = fpfr je Let (r ).pl(r vpt{r ) , Voir }d’r1, 


(IL.19) 


La méthode GGA donne des résultats en bon accord avec les résultats expérimentaux, elle est 
plus précise que la LDA avec une erreur réduite d’un facteur 3 à 5 comparés à la LDA[13] 
Effectivement, cette méthode est très efficace pour traiter les métaux de transition,[14] les 
réactions chimiques, [15,16] les solides et les systèmes magnétiques. [17,18] 

Il existe de nombreuses versions de la GGA telles que : Perdew, Wang et Becke [19] Perdew, 
Burk et Ernzerhof [20] (PBE). Wu-Cohen [21] ont développé la méthode WC-GGA, qui 
utilise sur une approximation simple et plus précise de l'énergie d'échange-corrélation, basée 
sur la fonctionnelle de corrélation de Perdew Burk Erzerhof (PBE). Bien que cette méthode a 
été conçue initialement pour les densités électroniques qui changent lentement, elle se révèle 


également efficace pour prédire les énergies d'échange des solides. 


ILS.3. L'approximation TB-mBJ: 


Les anomalies (sous estimations des résultats de calculs) observées sur les approximations 
LDA et GGA ont conduit Tran et Blaha à introduire une approche innovante en proposant leur 
propre potentiel d'échange, nommé BJ, [22] combiné avec le potentiel de corrélation de la 
LDA, pour avoir l’approche BJLDA. En se basant sur la BJILDA [23], Tran et Blaha [24] ont 
notablement amélioré leurs prévisions en modifiant le potentiel d'échange BJ pour créer une 
nouvelle méthode dénomée mBJ (Modified Becke Johnson). L'une des caractéristiques 
principales de la méthode TB-mBJ est sa capacité à fournir des valeurs précises du gap 
énergétique dans les semi-conducteurs. Cette adaptation relativement simple du potentiel BJ 
initial a permis d'obtenir une excellente concordance avec d'autres approximations de coûts 
plus élevés en raison de leur grande auto-cohérence, on cite à titre d’exemples les 
fonctionnelles hybrides [25,26] et la méthode GW [27, 28, 29]. L’expression du potentiel 
modifié BJ (mBJ) [30] proposé par Tran et Blaha en 2009 [24] se présente sous la forme 


suivante : 


Bj BR 1 5 2t (r) 
ve M=cv,, D+Gc-2) + = / @ (11.20) 
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m 
» 


Le symbole V: (r) fait référence au potentiel de Becker-Roussel (BR) [31], qui modélise le 
potentiel coulombien généré par les trous d'échange. Le terme pr) représente la densité 
électronique dépendante du spin, tandis que t_ (r) , Indique la densité d'énergie cinétique, le 


paramètre c (équation 11.20) est également pris en compte dans cette formulation. 


1 
2 
- 1_ pl) 73 
c= a + P 7. —J . é” (IL.21) 


Les paramètres a et y sont ajustables, leurs valeurs par défaut sont : à = - 0.012 et y = 1.023 


dd 
Bohr” , tandis que V ceu “Présente le volume de la cellule unitaire. 


IL. 6. La méthode des ondes planes augmentées et linéarisées 
(FP-LAPW) à potentiel plein: 


La méthode FP-LAPW, essentielle en physique des matériaux, est basée sur la DFT et utilise 
les ondes planes augmentées et linéarisées. Elle a été introduite dans les années 1930 et 
n’arrête pas d’être développée, elle a impacté positivement l’avancement dans le domaine de 
l’étude des propriétés électroniques des solides.[32] Elle combine la flexibilité des ondes 
planes avec la précision des ondes augmentées à l'intérieur des sphères atomiques, tout en 
linéarisant le potentiel électrique pour obtenir une description précise des interactions 


électroniques à l'échelle atomique.[33] 
IL.7. Méthode des ondes planes augmentées APW: 


La méthode APW (Augmented Plane Waves : APW), a été introduite initialement par Slater 
[34] en 1937, elle implique l'expansion des états propres d'une équation de Schrôdinger (11.2) 
pour une particule indépendante en utilisant des fonctions de base. Chaque fonction de base 
est représentée de manière unique dans les deux régions distinctes illustrées dans la Figure. 


112. 
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La région des sphères Muffin-Tin (MT) est une zone où les sphères sont disposées de telle 
manière à ne pas chevaucher et sont centrées sur les sites atomiques. Ces sphères notées 
"MT", ont un rayon R et à leur intérieur le potentiel a une symétrie sphérique et les fonctions 
d'onde sont radiales. 

Dans la zone interstitielle entre les atomes, où le potentiel est régulier, la fonction d'onde est 
représentée de manière à prendre en compte les variations graduelles qui relient les régions 
similaires à des atomes. 

Le potentiel est défini par la relation suivante : 


{V(r) = V(r), pour r > R} 
{V(r)=cst, pour r < R} ; (IL. 22) 


V(r) : est le potentiel Muffin-tin 


En conséquence, les fonctions de base pour les deux régions s’écrivent comme suit : 


{ fr) LA Y,, (r )- U,(r ,ÆE)}, pourr<R 


M) rs pourr>R (11.23) 
Q K 


{ dr) + > 2 «8 


Dr) Représente la fonction d'onde, Q désigne le volume de la maille unitaire, Fe sont les 


harmoniques sphériques, A Sont les coefficients associés au développement de la fonction 


l 


d’onde et U ; (r LE ,) st la solution radiale de l'équation de Schrôdinger (11.1) 
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Zone Interstitielle 
V=cCst 


Figure II. 2. L’espace associé à la maille unitaire du matériau simulé est divisé en deux régions : 


sphères autour des atomes et vide (espace interstitiel) entre les sphères. 


IL. 8. Méthode des Ondes Planes Augmentées Linéarisées 
(LAPW): 


La méthode des ondes planes augmentées linéarisées (LAPW) est une méthode de calcul 


développée par O. K. Andersen en 1975 [35]. Elle consiste premièrement à construire une 
fonction d'onde élémentaire U (r ,E D pour une énergie donnée FE . ensuite on effectue le 


développement de Taylor. Cette technique mathématique permet d'approximer une fonction 


inconnue à partir de ses dérivées : 


U(E,r) =U(E,r)+(E-E, )U'(E,r) +0 ((E : E)), (11.24) 
Les fonctions de base à l'intérieur des sphères sont exprimées sous forme d'une combinaison 
linéaire des fonctions radiales U : (r_,E : ) et de leurs dérivées par rapport à l'énergie U ; (r” 


JE d) , comme indiqué dans la relation suivante : 


fr) ZA .U, ( )+B 0 M]Y, @), pourr<R } 


Im 
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{D{r”) De Dee LE) SrSRt. 25 
G 


Les coefficients B " correspondant aux fonctions Ur) sont introduits. Pour déterminer À " 
et B im * On doit assurer la continuité à la surface des sphères Muffin tin. La dérivée de la 


fonction radiale U-/(r) est une expression générale obtenue par l’ensemble d'énergie E, de 


chaque atome avec une seule valeur fixe de E,. 


{- mn +440 + vo-E) rU GEJU (CE), (11.26) 


Il convient de noter que par rapport à la méthode APW, Le problème d'asymptote est 


définitivement éliminé, car même si U L SSt égal à zéro à la surface de la sphère, U | nese pas 


nulle . 

Il convient de mettre en avant que les électrons de cœur, totalement confinés au sein de la 
sphère MT, sont définis comme des états fondamentaux. D'un autre côté, les électrons restants 
sont qualifiés d'électrons de valence, existant à l'extérieur de la sphère MT tout en s'associant 
à d'autres atomes. Néanmoins, pour de nombreux éléments, établir une distinction nette entre 
ces catégories d'électrons s'avère complexe. Certains états ne relèvent ni des états 
fondamentaux n1 des états de valence et sont ainsi caractérisés en tant qu'états intermédiaires. 
Ces états partagent le même nombre quantique angulaire / que les états de valence, en dépit de 
leur nombre quantique principal n qui soit inférieur. Par conséquent, l'utilisation d'une énergie 
pour identifier les états partageant la même valeur de / devient une problématique. Ce 


dilemme trouve sa résolution avec l'introduction des orbitales. 


IL. 9. La méthode APW + lo: 


La méthode (APW+lo), intégré les fonctions d'onde locales (lo) pour une description 
améliorée des électrons à proximité des noyaux atomiques [35]. Sjôstedt et al. Ont démontré 
que la méthode de linéarisation de LAPW n'était pas optimale et pouvait être perfectionnée en 
exploitant quasiment toutes les fonctions de base dans les sphères, de manière similaire à la 


méthode APW originale (qui ne comportait qu'une fonction radiale, u). Dans cette variante 
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APW, l'énergie est fixée de la même manière que dans LAPW, utilisant une énergie 


prédéterminée E [36]. 


Les fonctions de base des ondes planes augmentées et des orbitales locales sont exprimées 


comme suit : 


{(r)= 24. U(r,E,.) +B 0 (GE, T0} pourr<R 


Im 


{D(r 7) = 0}, pour r >R (11.27) 
Pour introduire davantage de flexibilité dans les fonctions radiales, une nouvelle orbitale 
locale (lo) [36], distinguée par des lettres minuscules pour la différencier de la méthode 


LAPW, est ajoutée : 


of} EURE) +BU Ge E)1Y, @. (1.28) 


Il a été démontré que cette approche alternative converge rapidement vers les mêmes résultats 


que la méthode LAPW[37] 


IL. 10. Méthode LAPW+LO: 


La méthode LAPW+LO, constitue une extension de la méthode APW augmentées linéariser 
(LAPW), elle introduit les fonctions d'onde locales (LO) pour une amélioration significative 
de la description des états électroniques près des noyaux atomiques. Cette approche combine 
la flexibilité des ondes planes augmentées avec la précision des orbitales locales pour 
connaitre de manière plus précise le comportement électronique à proximité des noyaux 


atomiques. [38,39,35] 


La méthode LAPW+LO repose sur le principe fondamental des ondes planes augmentées 
(APW), où les fonctions d'onde électroniques sont développées en termes d'ondes planes à 
l'intérieur de chaque sphère atomique, tout en utilisant des fonctions radiales pour décrire les 


variations locales : 


tr JE >|4, Ur. E,)] c pt (e E, 1) ee % G E,,) Fe @)}, pourr <R 
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{fr ) = 0} pourr>R (IL.29) 


Les fonctions d'onde locales (LO) permettent d'affiner cette description en tenant compte de 


l'influence locale plus complexe des noyaux atomiques sur les électrons environnants. 


IL.11. Code de calcul Wien2k: 


Le code de calcul wien2k [40 ,41] est une série de programmes informatiques, utilisé pour le 
calcul des structures électroniques, 1l est basé sur la DFT. Plus précisément, il est conçu pour 
effectuer des calculs de type LAPW (Linear Augmented Plane Wave). Par la suite, ces 


informations peuvent être utilisées pour déduire diverses propriétés des matériaux, notamment 


+ Les positions atomiques après l'optimisation de la structure. 
* Les transitions de phase. 

+ Les propriétés électroniques 

+ Le magnétisme 


+ Les propriétés optiques, élastiques, etc... 


Les programmes principaux (ou "noyau") dans le code wien2k, essentiels pour chaque calcul, 
figurent sur le schéma illustré dans la Figure II. 3. Ces programmes sont classés en deux 
groupes distincts : 

Initialisation : 
NN, SGROUP, SYMMETRY, LSTART, KGEN, DSTART 

Cycle SCF : 
LAPWO, LAPWI, LAPW2, LCORE, MIXER 
Au cours de la phase d'initialisation, la structure d'entrée (fournie via le fichier .struct ) est 
traitée par les programmes NN, SGROUP, et SYMMETRY. NN, elle calcule les distances 
jusqu'aux voisins suivants et vérifie l'absence de chevauchements des sphères. 
SGROUP détermine le groupe spatial et les groupes ponctuels des sites d'atomes non 
équivalents dans le matériau. 


SYMMETRY génère des matrices pour les opérations de symétrie. 
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SCF CYCLE INITIALIZATION 
Pin 


LAPWO Input: Structure 
. lattice parameters 
enerates potential d 
8 p atoms(types, 


positions, Ruyr) 


V=Vc+VXc 


LAPW1 
solves eigenvalue 
problem, calculates 

valence bands 


checks for 
overlapping atomic 
spheres 


LCORE 
computes core states 
and densities 


SGROUP 


finds space group 
point groups of 
atoms 


LAPW2 
computes valence 
densities 


SYMMETRY 
generates symmetry 
operations 


Pval 


Pcore, Ecor 


MIXER creates 
new Pin by 
mixing (Pval+Pcore 
}and 


LSTART 
generates electronic 
densities of "free" 


Pold Pold atoms 


Converged ? generates k-mesh 


Pout 


DSTART 
superposes atomic 
densities 


Figure IL. 3. Organigramme d'un calcul général Win2k comprenant l'initialisation et Cycle SCF 


IL. 12. Conclusion: 


Le Chapitre II, dédié aux méthodes de calculs théoriques, contient des parties qui montrent les 
formalismes et outils mathématiques nécessaires à la compréhension des systèmes quantiques 
complexes. Beaucoup de concepts sont décrits tels que l'équation de Schrôüdinger, les 
méthodes développées pour solutionner l’équation de Schrodinger, la théorie de fonctionnelle 


de la densité (DFT) et la Méthode des Ondes Planes Augmentées et Linéarisées (FP-LAPW). 
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Ce chapitre évoque les avantages et les vertus des nouvelles méthodes de calculs développées 


sur l’avancée des travaux de simulation, appliquées à une large gamme de matériaux. 
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Chapitre LIT 


Etude théorique des propriétés structurales, électroniques, 
élastiques, optiques et magnétiques des composés Mg,,Cr,S 
(x= 0, 0.125, 0.25, 0.50, 0.75) dans la structure rock salt 
ferromagnétique. 
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Introduction: 


Grâce aux propriétés exceptionnelles des DMS du groupe II-VI et à leurs applications dans 
l'industrie, sur lesquelles nous avons précédemment discuté, nous avons entrepris une étude 
théorique approfondie des propriétés structurales, électroniques, élastiques, optiques et 
magnétiques du semi-conducteur MgS dopé au chrome (Cr) dans la structure de type rock-salt 
ferromagnétique. Les composés traités ont la composition chimique suivante : Mg,,Cr,S (x= 
0, 0.125, 0.25, 0.50, 0.75), où x représente le pourcentage de l'impureté (Cr) dans le composé. 
Les propriétés du composé sont calculées via le code Wien2k, qui utilise la méthode 


FP-LAPW. 


III. 1. Détails de calculs : 


Les propriétés des composés Mg,.,Cr,S (x = 0,0, 0.125, 0.25, 0.50, 0.75) sont calculées dans 
la structure rock salt (B1) ferromagnétique est effectué à l'aide de la méthode des ondes 
planes augmentées linéarisées à potentiel complet avec les orbitales locales FP LAPW + lo 
qui est mise en œuvre dans le code WIEN2K [1] en conjonction avec la théorie de la 


fonctionnelle de la densité de spin polarisée SPDFT [2]. Nous utilisons la structure RS des 


composés M£sS avec le groupe spatial N°225 (Fm3m) où les atomes Mg et S sont situés aux 
positions (0, 0, 0) et (1/2,1/2,1/2), respectivement. La substitution des atomes de Mg par des 
atomes de Cr est effectuée en fonction de la composition souhaitée. Le calcul des propriétés 
des composés est effectué en choisissant 16 atomes créés à partir d'une supercellule 2X2X2 et 
d'une supercellule cubique 1X1X1 contenant 8 atomes (voir la figure. III. 1). Les propriétés 
structurales et élastiques sont calculées avec l'approximation du gradient généralisé proposée 
par Wu et Cohen (WC-GGA) [3] , quant aux propriétés électroniques, optiques et 
magnétiques elles sont calculées via la méthode TB-mBJ (Tran-Blaha modified 
Becke-Johnson potential) [4], combinée avec la méthode PBE-GGA (approximation du 


gradient généralisé développée par Perdew Burk Erzerhof) [5]. 
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Le paramètre R,;K... est égal à 8 et les valeurs R,,- des atomes de Meg, Cr et S sont 
respectivement égales à 2.3, 2.4 et 2.01 (a.u.). Les états électroniques (3s°2p°) de Mg, (4s° 
3d*) de Cr et (3s/3p°) de S sont considérés comme des états de valence. 

L'expansion de Fourier de la densité de charge G.,, = 16 Ry, L,,, = 10, l'intégration de la zone 
de Brillouin dans l'espace K est faite avec 10X10X10 points pour les propriétés calculées et 
3000 points- k pour les propriétés optiques, le calcul auto-cohérent SCF est atteint à 10* Ry 


pour l'énergie totale et à 0.001 e pour la densité de charge des électrons. 
III. 2. Tests comparatifs d'approximation : 


Dans la méthode FP LAPW + lo, plusieurs approximations du potentiel d’échange-corrélation 
sont testées pour calculer le paramètre du réseau du composé binaire MgS. Le Tableau IIL.1 
présente différents tests effectués via plusieurs méthodes de calcul. 

Ces tests visent à trouver la valeur de "a" la plus proche possible de la valeur expérimentale a 
= 5.19 À. Comme illustré dans le Tableau IIL.1, l'approximation du gradient généralisé 


proposée par Wu et Cohen (WC-GGA) est celle qui répond le mieux à notre objectif. 


Tableau IIL. 1 Paramètre de réseau a (À) calculé avec différentes approximations. 


Approximation PBE-GGA | PBEsol-GGA | WC-GGA | Valeur Expérimentale 
66 


a(À) du MgS | 5.142 5.240 5.192 5.198 


Certaines approximations donnent des valeurs proches de la valeur expérimentale telles que 
les approximations WC-GGA et PBEsol-GGA ; la LDA, par contre, donne une valeur qui 


présente un écart signifiant par rapport à la valeur expérimentale. 


66 


IL. 3. Propriétés structurales : 


Les propriétés structurales des composés Mg,.,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75) dans la 
phase rock salt (RS) ferromagnétique reposent sur le calcul du paramètre de réseau à 
l'équilibre « a, », le module de compressibilité B et sa première dérivée B'. Les valeurs de ces 
paramètres, affichées dans le tableau. IIT. 2, sont estimées en utilisant l'approximation 
WC-GGA par ajustement des courbes donnant l'énergie totale en fonction du volume à 
l'équation d'état de Birch-Murnaghan [6]. 

Pour le M£gS pur, les valeurs calculées de a,, B, et B’ sont en bon accord avec les valeurs 
expérimentales et sont similaires aux valeurs théoriques de la littérature. Pour les composés 


ternaires, les valeurs des paramètres structuraux ne sont pas disponibles dans la littérature. 
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Figure IL. 1. Structures cristallines des composés (a) Mg s75Cro 1258, (b) Mgo5Cro25S, (C) Mg 50CTo 509 


et (d) Mgo25Cro758. 


Tableau IIL. 2 Valeurs calculées du paramètre de réseau a (À), du module de compressibilité B (GPa), 
et de sa dérivée première B', des énergies minimales E;y(ev), Earm (eV), Exu(ev) dans les trois phases 
FM, AFM et NM, et de l'énergie de formation E; (ev) des composés Mg,,Cr,S (x = 0,0, 0.125, 0.25, 
0.50 et 0.75). Toutes les valeurs sont obtenues via l'approximation WC-GGA et comparées à d'autres 


valeurs théoriques et expérimentales de la littérature. 
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Composé Notre travail Valeur Autres valeurs B(GPa) Eru(ev) Exrm(eV) Exu(eV) E{eV) 
a(A°) Expérimentale théoriques 
a(A°) a(A!) 


À 5.189', 79.21 4.1 
5.18°,5.16", | 81077.712 ù 
Structure RS 5.142! 


Mgo.s7sCro125S 
MgorsCro2sS 
Mgo.s0Cro.50S 


Mgo2sCrorsS 


On constate que l'élévation du taux du chrome (Cr) dans les composés entraîne une 


diminution de la valeur du paramètre de réseau a, et une augmentation de la valeur du module 
de compressibilité B indiquant une tendance croissante de la dureté des composés. Les valeurs 
de l'énergie minimale des composés Mg,,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75) dans la 
structure RS a été calculée dans les trois phases : non-magnétique (NM), ferromagnétique 
(FM) et antiferromagnétique (AFM) sont présentées dans le tableau. IIL. 2. Pour chaque 
composé, la valeur la plus faible de cette énergie est celle de la phase ferromagnétique, 
indiquant que les composés sont stables dans cette phase. En résumé, ces résultats confirment 
la stabilité des composés ternaires dans la structure rock salt ferromagnétique. L'énergie de 
formation des ternaires a également été calculée, les valeurs sont indiquées dans le tableau. 
II. 2, elles sont toutes négatives, ce qui indique la possibilité de leur élaboration et confirment 
également leur stabilité structurale. Le composé le moins dopé (x = 0.125) a la valeur la plus 
négative, ce qui en fait de lui le meilleur composé en termes de stabilité. Rappelons que 
l’énergie de formation des composés a été calculées depuis la formule suivante [13] : 


E;= Eugi-xCxs- (n,E"8 + n,ET+n,E°), 
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Où E"* est l'énergie totale du composé contenant x% de Cr, Eu, Et et Es sont les énergies 
totales des éléments respectifs Mg, Cr, et S et n, n,et n, les nombres d'atomes respectifs des 


éléments Mg, de Cr et de S dans les supercellules. 


IL. 4. Propriétés électroniques : 
IL. 4. 1. Structures de bandes : 


Les bandes d'énergie des composés Mg,,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75), calculées par 
la méthode TB-mBJ dans la phase RS ferromagnétique, sont présentées dans les figures IT. 2. 
à 6. La Figure IT. 2 montre une bande interdite indirecte (1'-X) pour le M£S pur, les quatre 
composés ternaires ont des bandes interdites directes dans la direction I de haute symétrie 
dans la première zone de Brillouin et dont la valeur diminue avec l'élévation du taux de Cr 
dans les composés. On remarque que le haut de la bande de valence reste stationnaire tandis 
que le bas de la bande de conduction se déplace considérablement vers le niveau de Fermi 


quand le taux de Cr augmente. 
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Figure IIL. 2. Bandes d'énergie de MgsS pur, calculées par l'approximation TB- mBJ. 
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Figure IIL 3. Bandes d'énergie polarisées en spin du ternaire Mg s5Cro5S, calculées par 


l'approximation TB- mBJ. 
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Figure LIL. 4. Bandes d'énergie polarisées en spin du ternaire Mg 5Cro2sS, calculées par 
l'approximation TB- mBJ 
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Figure IL. 5. Bandes d'énergie polarisées en spin du ternaire Mg, s0Cr0.508 (a), calculées par 


l'approximation TB- mBJ. 
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Figure IIL. 6. Bandes d'énergie polarisées en spin du ternaire Mg, ,.Cro-sS, calculées par 


l'approximation TB- mBJ. 


Le résultat le plus remarquable est le changement de la nature du gap du composé binaire 


Me&s par l'effet de son dopage au chrome (Cr) où, en effet, on passe d'un gap indirect dans le 


M&sS ![14,12,10,15] à un gap direct dans les composés Mg,.,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 


0.75). Dans les quatre composés ternaires, on observe la présence d'états électroniques au 


niveau de Fermi dans le canal des spins majoritaires indiquant le caractère métallique; par 


contre, dans le canal des spins minoritaires, l'existence d'une bande interdite confirme le 


comportement semi-conducteur. 


La polarisation en spins au niveau de Fermi (P) est un bon indicateur de la demi-métallicité 


des (DMS), son expression est la suivante : 


. NT(E,)-M (E,) 


MATE) 


(IL. 1) 
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Où NT(E,) et N(E,) sont respectivement, les densités de spin-up et de spin-down au niveau de 
Fermi. La valeur de la polarisation en spin au niveau de Fermi est maximale (100 %) dans 
tous les composés ternaires. 

Les valeurs calculées de la bande interdite E, et du gap demi-métallique Gi des différents 
composés sont indiquées dans le tableau IL. 3. Gyn, représente le minimum entre l'énergie la 
plus basse Ec,;,(T) du spin majoritaire ou minoritaire dans les bandes de conduction et les 
valeurs absolues de l'énergie la plus élevée du spin majoritaire ou minoritaire dans la bande de 
valence Ev,,,, (1). Dans notre cas, la valeur de Gyn, est obtenue à partir de la différence entre le 
minimum de la bande de conduction et le niveau de Fermi dans le canal des spins 


minoritaires. 


Tableau IIL. 3. Valeurs calculées de (Eg), de (G:m) et des moments magnétiques totaux et locaux des 


composés Me..,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75) par la méthode TB-mBJ. 


js _. nue 


M£os75Cro.5S 
Mgo.75Cro25S 


Concernant le Mes pur, la valeur théorique de E, relevée de la littérature est (4.010 eV) [14] 


est proches de la nôtre (4.043 eV), l’écart entre les deux valeurs est faible et vaut 0.033 eV. 
Dans les composés respectifs Mgss75Cro125S, MgosCtro2sS, Mgos0Cto50S et Mgo5Cto75S, les 
valeurs de Eg sont égales à 3.86 eV, 3.41 eV, 2.73 eV et 2.65 eV, on remarque une 


décroissance de la valeur de Eg en fonction du taux de chrome. 
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Les valeurs du gap demi-métallique Gym (tableau IIL.3) diminuent aussi avec le taux de Cr 
dans les composés mais jusqu’à la concentration Cr-50 % ; en effet, dans le composé le plus 
riche en chrome (75%) la valeur du gap subit une augmentation. Rappelons que la valeur 
maximale de Gr est celle du composé à faible teneur en chrome (12,5% Cr). 

Le tableau IIL. 3 montre que le moment magnétique total de tous les composés ternaires est un 
nombre entier et vaut 4 13, il est réparti entre les atomes qui forment le composé et il est 
important pour le cas du chrome (Cr) comparativement avec les moments des atomes non 


magnétiques Mg et S et de l'espace interstitiel. 


III. 4. 2. Densité des états électroniques: 


Les figures IIL. 7 à IIL11 montrent les courbes de densité d'états totale TDOS et partielle 
PDOS des composés Mg,,Ci,S (x = 0.0, 0.125, 0.25, 0.50, 0.75) dans la structure rock salt 
ferromagnétiques. Les courbes TDOS et PDOS du Mes (Figure IIL. 7.) montrent une forte 
densité des états de type p-S au sommet de la bande de valence près du niveau de Fermi, la 
bande de conduction est, par contre, occupée par les états s-Mg, p-Me, s-$, et p-S. Les figures 
IT. 8. (a), 9. (a), 10. (a), 11. (a) indiquent le comportement demi-métallique des composés 
ternaires, leur polarisation de spin au niveau de Fermi est de 100%. L'hybridation p-d 
impliquant les états 3p-S et 3d-Cr est observée près du niveau de Fermi, elle est d'autant plus 
importante que le pourcentage de Cr dans les composés est élevé. 

Dans les deux canaux de spins majoritaires et minoritaires, les courbes de densité d'état des 


composés sont asymétriques en raison des fortes hybridations p-d suscitées. 
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Figure IL. 7. Courbes de densité d'états totale et partielles du composé binaire Mgs, calculées à l'aide 
de l'approximation TB-mBJ. 
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Figure IIL. 8. Courbe de densité d'états totale (a) et partielles (b) (c) du composé Mg, 475Cr0 125, 
calculées à l'aide de l'approximation TB-mBJ. 
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Figure IIL. 9. Courbes de densité d'états totale (a) et partielles (b) (c) du composé Mgs75Cro25S, 


calculées à l'aide de l'approximation TB-mBJ. 
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Figure II. 10. Courbes de densité d'états totale (a) et partielles (b) (c) du composé Mg, s0Cro508, 
calculées à l'aide de l'approximation TB-mBJ. 
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Figure IIL. 11. Courbes de densité d'états totale (a) et partielles (b) (c) du composé Mgs25Cro7sS, 
calculées à l'aide de l'approximation TB-mBJ. 


Les courbes TDOS et PDOS (Fig. IL. 8. b-c, Fig. 9. b-c, Fig. 10. b-c, Fig. 11. b-c), montrent 
que les sommets des bandes de valence des composés sont formés par les états électroniques 
3d-Cr et par une faible contribution des états 3p-S et 3p-Mg au bas des bandes de valence 
situées entre -6.93 eV et 2.5 eV pour MgssCronsS, entre -7.02 eV et 2.5 eV pour 
Mgo5Cro2sS, entre -7.39 eV et -2.57 eV pour Mg, :0Cr059 et entre -6.99 eV et -1.41 eV pour 
Mgo2sCto7sS. 

Dans le canal des spins minoritaires, la bande de conduction contient les états 3d-Cr ainsi que 


les états p-Mg et p-S mais avec une faible représentation. Sous l’effet de la symétrie 
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octaédrique du champ cristallin et par suite aux interactions d'échange entre les orbitales 
3d-Cr et p-S, l'orbitale 3d-Cr est divisée en deux : les orbitales localisées triplement 
dégénérées t,, (dyz, dxz et dxy) de basse énergie et les orbitales doublement dégénérées de 
haute énergie e, (dx°-y”, dz”). La figure IIL.12 présente le champ cristallin octaédrique, où six 
ligands entourent l'ion métallique, les orbitales d se divisent en deux ensembles de niveaux 
d'énergie différents. Les trois orbitales les moins énergétiques (dyz, dxz et dxy) sont appelées 
collectivement l'ensemble "t,,", tandis que les deux orbitales les plus énergétiques (dx°-y° et 


dz”) forment l'ensemble "e,". 


Energie x°-y° ee, 
1° r 


er 


" 
e 


Figure III. 12. La dégénérescence des orbitales « d » en orbitales t,, et e, sous l'action du champ 


cristallin octaédrique. 
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IL. 5. Propriétés élastiques : 


Compte tenu de la symétrie cubique des composés Mg,,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 
0.75), trois constantes élastiques (C1, C1, Cx) suffisent pour décrire leurs propriétés 
élastiques. La méthode de calcul utilisée pour l’obtention des trois constants élastiques est la 
méthode FP-LAPW [16]. 

Le tableau IIT. 4 regroupe les valeurs, à pression nulle, des constantes élastiques C;,, C;, et 
C4, elles sont comparées à d'autres valeurs théoriques antérieurement publiées [11,17]. Pour 
le Mes pur, les valeurs des constantes C.;, et C;, sont relativement plus élevées, l'inverse est 
observé pour la valeur de C;, ; les différences observées sont attribuables à la méthode de 
calcul GGA [18]. 


Tableau LIL. 4. Valeurs des constantes élastiques des composés M£s et Mg, ,Cr,S (x = 0.125, 0.25, 


0.50, 0.75), calculées en utilisant l'approximation WC-GGA. Les valeurs obtenues sont comparées aux 


valeurs théoriques de la littérature. 
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La stabilité mécanique des matériaux cubiques doit répondre aux critères de stabilité de 
Born-Huang [19,20] qui sont les suivantes : 


Ca > 0 
(Cun-Cr2)>0 
(Cu+2C12)>0 


Dans notre cas, les valeurs C;,, C;, et C;, des composés Me,.,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 
0.75), vérifient les conditions de stabilité de Born-Huang ; par conséquent, les composés sont 
considérés élastiquement stables. Les valeurs de C;, étant supérieures à celles de C;, et C4, 
les composés ont donc une résistance à la compression unidirectionnelle [21]. 

L'anisotropie élastique À des composés est estimée à partir de l'expression (2) ci-dessous, elle 
est égale à l'unité (A=1) pour les composés isotropes et elle est soit supérieure soit inférieure à 
l'unité pour les composés anisotropes [22]. Les valeurs de A affichées dans le tableau II. 5 
indiquent que le MgS est isotrope (A = 1) et que les composés ternaires Mg,,Cr,S (x = 0.0, 
0.125, 0.25, 0.50, 0.75) sont anisotropes étant donné que leurs valeurs de A sont différentes de 
l’unité. On remarque que l’anisotropie des composés est fortement dépendante de la 
concentration de l’impureté (Cr) introduite dans la matrice du Mgs. Les valeurs de grandeurs 
Cs, Gun, et E (équations IIL 3, 5 et 8 respectivement) subissent une diminution avec 
l’augmentation du taux de Cr, celles des grandeurs B/G; et v (équations IIL. 4, 5 et 9 
respectivement), par contre, augmentent. 

La ductilité (fragilité) du matériau est évaluée sur la base du rapport Pugh B\G; [23]. Un 
matériau est considéré comme ductile si B/G;; est supérieur à 1.75 et comme fragile si B/Gx 
est inférieur à 1.75. 

Les valeurs de B/Gx; dans le tableau II. 5. indiquent que les composés ayant des 
concentrations x = 0.0, 0.125, 0.25 et 0.50 sont fragiles et que le composé contenant le 


maximum d’impureté (x = 0.75) est ductile. 
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Les valeurs de ce paramètre sont indiquées dans le tableau II. 5. 


2C,, 
A = rer (HI. 2) 
C..—C 
C. - ( ne é] (IL. 3) 
C..+2C 
B = se) - = u) (HI. 4) 
G +G 
G, = = : (LIL. 5) 
C.—C .+3C 
G,= —_— — (III. 6) 


CR 4C,,+3(11 —C,) QU. 7) 


9BG 
E= (II. 8) 
__ 3B-2G 
Ve Se (HI. 9) 


Tableau IIL 5.Valeurs calculées du facteur d’anisotropie élastique A, du module d'onde de 
cisaillement Cs, du module de cisaillement G, (GPa), du module de Young E (GPa), du coefficient de 


Poisson v et du rapport B/G;, en utilisant l'approximation WC-GGA. 
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MgorsCro2s$ 62.95 51.24 125.23 | 0.22 1.46 
Mgo.50C 0.508 52.21 42.78 107.16 0.25 | 1.69 


MgozsCro7sS 47.15 34.93 88.78 0.27 | 1.85 


Les figures III. 12, 13, 14, 15 et 16 représentent les formes géométriques à 2D et à 3D des 
modules B, E et G;; des composés Mg,.,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75). A 2D, elles 
sont représentées dans les plans (100), (010), (001). Les figures de B ont la forme sphérique 
synonyme d'anisotropie, celles de E et G;; ont une forme différente de la sphère indiquant leur 


anisotropie. 
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Figures IIL. 12. Représentation 3D de la dépendance directionnelle (a) du module de compressibilité, 
(b) du module d'Young et (c) du module de cisaillement du composé binaire 
Mgs. 
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Figures LIL. 13. Représentation 3D de la dépendance directionnelle (a) du module de compressibilité, 
(b) du module d'Young et (c) du module de cisaillement du composé ternaire Mg 875Cro 1258. 
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Figures III. 14.Représentation 3D de la dépendance directionnelle (a) du module de compressibilité, 
(b) du module d'Young et (c) du module de cisaillement du composé ternaire Mg, ;5Cro2sS. 
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Figures IIL. 15. Représentation 3D de la dépendance directionnelle (a) du module de compressibilité, 
(b) du module d'Young et (c) du module de cisaillement du composé ternaire Mg, :5Cro508. 
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Figures IIL. 16. Représentation 3D de la dépendance directionnelle (a) du module de compressibilité, 
(b) du module d'Young et (c) du module de cisaillement du composé ternaire Mg, ,Cro7sS. 
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III. 6. Propriétés optiques : 


La constante diélectrique complexe e(æ) donne des informations sur les transitions 
électroniques possibles dans un matériau lorsqu'il interagit avec un rayonnement 
électromagnétique, elle s’écrit comme suit : 

e(@) = &, (©) He (), où i= (-1)'? est le nombre complexe. 

La partie réelle £,(w) est liée à la polarisation du milieu et la partie imaginaire £,() est en 
étroite relation avec la réponse du matériau à l‘excitation de l’onde électromagnétique, les 
deux grandeurs sont liées par les relations de Kramers-Kroning [24]. Les calculs des 
grandeurs optiques sont effectués via l’approximation TB-mBJ dans l’intervalle d’énergie 
entre 0 et 40 eV 

La partie réelle de la constante diélectrique &,(@) (Fig. IL. 17. a) informe sur la métallicité des 
composés et sur leur pouvoir à réfléchir la lumière incidente. Au-delà de 13 eV, les courbes 
e(o) relatifs aux composés Mg,,Cr,S (x = 0.00, 0.125, 0.25, 0.50, 0.75) ont des allures 
similaires et des amplitudes presque égales, on n’observe pratiquement par l’effet du dopant, 
les courbes sont superposées. L’effet du dopant est, par contre, observé sur le reste du spectre 
optique #,() des composés notamment dans la région de la lumière ultraviolette au-delà de 6 
eV où on note une baisse considérable de l’amplitude de £,(&) et un profil totalement différent 
de celui observé aux grandes énergies (E>13 eV). Dans l’intervalle énergétique entre 13 eV et 
16 eV, £&,(w) est négatif indiquant le comportement métallique et la haute réflectance des 
composés, au-delà de 16 eV £,(@) est positif et son amplitude est faible. 

Aux basses énergies, entre 0 et 12 eV, l’amplitude de £,() est considérable. Un pic intense est 
localisé autour de 6 eV dans le M£S pur et dans les composés, son intensité et sa position 


énergétique sont dépendantes du taux de chrome. Dans cette région, le signe de £,(w) des 
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composés est positif entre 0 et 6 eV, il s’annule dans chaque composé à une énergie située 
dans l’intervalle de 7 eV à 9 eV et au-delà de laquelle il devient négatif, il s’annule une 
deuxième fois plus loin aux entours de 17 eV et devient positif vers la fin du spectre, au-delà 
de 17 eV son amplitude varie d’une manière monotone avec l’énergie et elle est insensible au 
pourcentage du dopant dans les composés. 

Les valeurs statiques (£:(0)) des composés respectifs MgS, Mgo5Cto15S, Mgo5Cto258, 
Mg 0Ct0.508 et Mgo5Cr9758 sont les suivantes : 4.64, 2.95, 3.03, 3.22 et 3.24. Rappelons que 
pour les composés ternaires les valeurs de £,(0) sont prises dans le canal des spins 
minoritaires (spin-down) là où le composé possède le caractère semi-conducteur. Remarquons 
qu'avec l’augmentation de la concentration de Cr, la valeur de £,(0) des semi-conducteurs 
ternaires augmente et inversement la valeur de leurs gaps énergétiques diminue. Dans notre 
cas, les valeurs de £,(0) et du gap énergétique E, (voir l'équation III. 10) satisfont au modèle 
de Penn [25] selon lequel ces deux grandeurs sont en première approximation liées par 
l’équation (III. 9) ci-dessous. 

e(0) = 1+ (ho,/E,) (HI. 10) 

La figure IIL. 17. b représente la partie imaginaire de la constante diélectrique £,(©) du 
composé binaire M£gS pur et des composés ternaires Mg,.,Cr,S (x = 0.125, 0.25, 0.50, 0.75), 
les courbes sont tracées entre 0 à 40 eV et renseignent sur les transitions électroniques 
provoquées suite à l’absorption des ondes électromagnétiques par les composés. Lorsqu’on 
passe du MgS pur aux composés, les courbes £,(©) changent spectaculairement. Le 
changement fondamental est observé dans la région des énergies ultraviolettes (E<S eV) où 
on note la présence de pics supplémentaires quasiment absents dans le semi-conducteur 
binaire MgsS, les plus importants sont situés autour de 5 eV et en dessous de 1 eV, un 
troisième pic de faible amplitude est présent autour de 3 eV. Les pics caractéristiques du 


composé binaire MgS pur sont situés dans la région entre 6 eV et 10 eV, et il est facile de voir 


93 


qu'ils sont affectés par le dopage. En effet, l’augmentation du taux de Cr conduit à un sensible 
déplacement des positions des pics vers les grandes énergies et à une baisse considérable de 
leurs intensités. L'origine de ces pics, si on se réfère aux courbes TDOS et PDOS de M£gsS pur, 
résulte des transitions électroniques mettant en jeu les états électroniques de type p-S occupés 
dans la bande de valence et les états électroniques s-Mg inoccupés dans la bande de 
conduction. Les pics supplémentaires sont dû au métal Cr, le pic entre 0 et 2.5 eV est dû aux 
transitions électroniques impliquant les états occupés 3d-eg/Cr présents dans le haut de la 
bande de valence des composés ternaires et les états électroniques vides p-S situés dans le bas 
de la bande de conduction dans le canal des spins majoritaires, le pic autour de 5 eV provient 
des transitions entre les états p-S occupés existant dans le haut de la bande de valence et les 
états 3d-Cr vides se trouvant dans le bas de la bande de conduction dans le canal des spins 
minoritaires. Avec l’augmentation du taux de Cr, les deux principaux pics caractérisant la 
région entre 0 et 6 eV gagnent en amplitude et leurs positions sont légèrement recalées vers 


les faibles énergies. 
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Figure IL. 17. Partie réelle #, (&) (a) et partie imaginaire £&, (©) (b) de la fonction diélectrique des 
composés Mg,,Cr,S (x =0.00, 0.125, 0.25, 0.50, 0.75), calculées à l'aide de l’approximation TB-mBJ 
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III. 6.1 Indice de réfraction n (o) et indice d'extinction k (@): 

L'indice de réfraction n(@) renseigne sur la propagation d’une onde électromagnétique dans 
un solide, sa connaissance est fondamentale pour la caractérisation des matériaux utilisés dans 
les dispositifs optoélectroniques [10]. La figure IIL 18. a représente la fonction n(w) dans 
l'intervalle entre 0 et 40 eV, son allure indique que les composés sont des milieux fortement 
dispersifs car n(®) varie considérablement en fonction de l’énergie de l’onde 
électromagnétique, en particulier entre 0 et 16 eV. L’effet du dopant sur la courbe n(œ) est 
bien marqué dans la région de la lumière ultraviolette et dans la région des basses énergies 
entre O et 2.5 eV. Dans la région de la lumière ultraviolette la dispersion de n(œ) est 
importante et on note un maximum dont la position est autour de 6 eV dans le M£sS pur et elle 
est légèrement au-dessus de 6 eV dans les composés ternaires, la valeur de ce maximum est à 
peu près égale à 3.5 dans le M£S pur et varie dans les composés selon la concentration du 
dopant. Au-delà de 6 eV n(w) décroit, sa valeur à 17 eV vaut approximativement 0.3. Pour les 
énergies tel que E>17 eV, n(®) remonte monotonement la pente et atteint la valeur 
approximative de 0.8 vers la fin du spectre optique. Dans la région du visible, la valeur de 
n(w) est autour de 2.5 aux environs de 3 eV. Au zéro des énergies, l’indice de réfraction 
statique prend les valeurs 2.15, 4.85, 3.00, 3.68, et 5.13 dans les composés respectifs Mgs, 
Mgos75Cto125S, Mgo7sCto25S, Mgos0CT0.508; Mgo25Cto75S, ces résultats montrent la sensibilité 


de l’indice statique au dopage par le chrome. 


L'indice d’extinction k(w) est une grandeur physique qui caractérise le pouvoir d’un matériau 
à freiner la propagation d’une onde électromagnétique. La figure IL. 18. b donne la variation 
de k(&@) en fonction de l’énergie de l’onde électromagnétique entre 0 et 40 eV. Les maximums 
de k(w) sont atteints dans la région des ultraviolets, deux pics intenses sont relevés dans le 
Mes pur autour de 7.5 eV et 10.0 eV d’amplitudes respectives égales à 2.6 et 2.4. Dans les 


composés ternaires, sous l’effet du chrome, les pics intenses subissent des changements, leurs 
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intensités diminuent et leurs positions se déplacent vers les grandes énergies. Aux basses 
énergies, les courbes k(w) des composés ternaires sont marquées par la présence de pics 
supplémentaires inexistants sur la courbe k(w) du composé binaire, leur présence dans cette 


région est liée étroitement à l’élément dopant Cr. 
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Figure III. 18. (a) Indice de réfraction n (œ) et (b) Coefficient d'extinction K (@) des composés 


Mg.,Cr,S (x=0.00, 0.125, 0.25, 0.50, 0.75), calculés via l'approximation TB-mBJ. 


IL. 6. 2 Conductivité optique : 

La figure IIL. 19. Montre les courbes de conductivité optique o(œ) = 2nc £,(w)/X des composés 
Me£.,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75), tracées toujours dans l’intervalle entre 0 et 40 eV. 
«À » est la longueur d’onde du rayonnement, «ec » est la vitesse de la lumière et £,(w) est la 
partie imaginaire de la constante diélectrique. Dans le semi-conducteur binaire MgS, comme 


dans les composés ternaires, les valeurs maximales de o(@) se situent dans la région des ondes 
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ultraviolettes où on relève deux pics de forte intensité positionnés aux énergies 
approximatives 7.5 eV et 9 eV dans le Mes pur. Les courbes 6(«) des ternaires présentent des 
écarts par rapport à celle de MgsS, ces écarts sont dus à l’ajout du chrome et concernent 
principalement les intensités et les positions des pics les plus intenses, on note une baisse de 
l'intensité des pics et un déplacement de leurs positions vers les hautes énergies. Aux basses 
énergies, 6(&) fait paraitre des pics de faibles intensités dont l’origine est associée aux états 
3d de Cr présents au voisinage du niveau de Fermi, leurs intensités et leurs positions varient 
avec la concentration du chrome. Au-delà de 10 eV, les courbes o(®) des composés 
décroissent avec l’augmentation de l’énergie des ondes électromagnétiques et leurs 


amplitudes ne présentent pas un écart signifiant. 
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Figure IIL. 19. Conductivité optique des composés Mg, ,Cr,S (x=0.00, 0.125, 0.25, 0.50, 0.75), 


calculées par la méthode TB-mBJ 


IL. 6.3 Coefficient d’absorption a(o) : 
Le coefficient d’absorption a(w) indique la capacité d’un matériau à absorber la lumière 


électromagnétique. Sur la figure. III. 20 sont représentées les courbes a(w) des composés 
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Meg.,Cr,S (x = 0.0, 0.125, 0.25, 0.50, 0.75) entre 0 et 40 eV. Il est clair que l’absorption des 
composés est importante dans le domaine des ultraviolets, l’absorption supplémentaire 
observée aux faibles énergies donne aux composés un profil a(w) qui diffère de ceux des 
composés ternaires Mg,,Cr,S (x = 0.125, 0.25, 0.50, 0.75). Des différences sont également 
reliées concernant l’aspect des courbes surtout aux basses énergies entre 0 et 6 eV. Entre 12 et 
20 eV, les courbes se rejoignent, une légère différence est observée au-delà de 20 eV où 


l’amplitude montre une sensibilité à l’augmentation du taux de Cr. 
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Figure IL. 20. Courbes a(w) des composés Mg,.Cr,S (x=0.00, 0.125, 0.25, 0.50, 0.75), calculées par 
la méthode TB-mBJ. 


III. 7. Propriétés magnétiques : 


Dans cette section, seront présentés les propriétés magnétiques des ternaires Mg,,Cr,S (x = 
0.125, 0.25, 0.50, 0.75). Nous donnons les valeurs calculées de leurs constantes d'échange 


Noa (s-d) et NB (p-d). 
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Les interactions d'échange entre le spin de l'impureté magnétique et le spin des porteurs de 
charge sont évalués à partir des constantes d'échange Na et Nop à l'aide de l'hamiltonien 
Kondo [26], et peuvent donc être calculées à partir de la répartition du spin entre les bandes 


de valence et de conduction comme suit : 


AE. 
N 0% © 

AE, 
A 7 x(s) 


Avec ÂË, — E. (1) = E. (T) et AE, = E, () = E, (D) 


x est la concentration de Cr et 
(s) est la magnétisation moyenne de Cr. 


Tableau IIL 6. Calcul de la division du bord de la bande de conduction (AË:), de la division du bord 
de la bande de valence (AË,) et des constantes d'échange (Noa et NB) de chaque site dans les 


composés Me... Cr,S (x = 0.125, 0.25, 0.50, 0.75), à l'aide de la méthode TB-mBJ. 


Composé AE, AE, Noa No 
(eV) (eV) 


Mgos75CTo125S -0.219 2.244 -0.483 4.937 


MgorsCro25sS | -0.727 1.467 -1.574 3.176 
M£o50Cro.50S | -1.612 2.036 -1.749 2.209 
Mpgo2sCro75S | -1.624 2.282 -1.156 1.625 


Le tableau IIL. 6. Montre que les valeurs des constantes d'échange ont des signes opposés où, 


en effet, N o% SSt négatif et N oP est positif. Cette situation s’explique par l'existence d’un 


couplage d'échange antiferromagnétique entre la bande de conduction et les états 3d du Cr, et 
d’un couplage d'échange ferromagnétique entre la bande de valence et les états 3d du Cr, 


respectivement. 
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Chapitre IV 
Propriétés optoélectroniques et mécaniques des 


semi-conducteurs MgXS (X= Cd, Zn) à large bande interdite 
dans la phase stable non magnétique. 


101 


Introduction : 

Dans cette partie, nous présentons une étude théorique spécifique des propriétés structurales, 
électroniques, mécaniques et thermoélectriques du composé binaire MgS dopé par le zinc 
(Zn) et par le cadmium (Cd). L’étude est effectuée en utilisant la méthode FP-LAPW, 
implémentée dans le code Wien2k et les composés sont étudiés dans la phase rock-salt (RS) 
non magnétique. Rappelons que les éléments Cd et Zn ont dans leurs couches externes une 
orbitale “d” stable et entièrement pleine (Cd-4d!° et Zn-3d!°), contribuant ainsi à leur stabilité 


chimique. 
IV. 1. Détails des calculs : 


Les calculs dans cette section sont effectués à l'aide du code Wien2k. Nous avons calculé les 
propriétés structurales et élastiques via la méthode WC-GGA, les propriétés électroniques et 
optiques via la méthode TB-mBJ et les propriétés thermoélectriques via le code Boltztrap, en 
fixant L,,,=10 et G,x = 16. 

Les fonctions d'onde dans l'espace interstitiel (espace inter-atomique) sont planes, avec un 
seuil fixé à Rur * Kmax — 8. Pour les atomes Cd/Zn, Mg et S, les valeurs de R, sont 
respectivement égales à 2.5, 2.3 et 2.04 (a.u). Les états (3s° 3p°) de Mg, (4s° 3d!°) de Zn, (5s 
4d'°) de Cd et (3s 3p°) de S sont considérés comme des états de valence. 

Le calcul a été effectué dans une supercellule cubique 2X2X2 composée de 16 atomes, le 
maillage choisi est de 1000 points-k pour toutes les propriétés étudiées sauf pour les 
propriétés optiques où on a utilisées 3000 points-k. Les critères de convergence auto 
consistants ont été appliqués avec une optimisation de l'énergie à 10* Ry et une densité de 


charge électronique à 0.001 Ry. 
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IV. 2. Propriétés structurales : 


Nous avons utilisé l'équation d'état de Birch-Murnaghan (BMES) [1] avec l'approximation 
WC-GGA pour étudier la structure cristalline la plus stable des composés M£s et 

Mgo875X0125S (X=Cd, Zn) en optimisant leurs courbes exprimant l'énergie totale en fonction 
du volume, les résultats sont présentés dans la figure IV.1. Les courbes indiquent une stabilité 


des composés ternaires MgossCdo2sS et MgossZnosS (Figures IV. 1.(b) et (c)) dans la 


structure rock salt (RS) de réseau cubique à faces centrées et de groupe d’espace Fm3m 
(225). Les structures cristallines des composés ternaires sont construites en remplaçant un 
atome de Mg par un atome de Cd (Zn) ; une supercellule primitive de 16 atomes est ainsi 
formée. Les atomes de Mg et de S occupent les positions respectives de Wyckoff (0,0,0) et 
(2, 2, V2), tandis que l'élément X= Cd, Zn occupe les sommets et les centres des faces du 
cube. Les valeurs de a (en À) et de B (en GPa) des composés sont présentées dans le tableau. 


IV. 1 et comparées aux valeurs précédemment publiées. 


La valeur de la constante de réseau du MgsS pur, égale à 5,19 À, concorde très bien avec la 
valeur expérimentale "a" = 5,19 À [2]. Celles des composés ternaires Mgog5CdosS et 
MgogrsZhosS sont respectivement égales à 5.23 À et 5.17 À. Il est également à noter que la 
valeur du module de compressibilité du composé Mg, s75Zno125S est supérieure à celles des 
composés MgsS et Mgos5CdosS. Les valeurs de l'énergie de formation du Mes, du 
Mgos75CdosS et du Mgos5Z55S figurant dans le tableau. IV. 1 sont toutes négatives, 
indiquant leur stabilité thermodynamique [3] et la possibilité de les synthétiser 


expérimentalement. Mgs75Zho0125S est le composé le plus stable, avec une énergie de 


formation (E -69.10 ev) plus négative que celle des deux autres composés. 


Tableau IV. 1. Paramètre de réseau a (À), module de compressibilité B (GPa), énergie de formation E; 
(eV) et énergie du Gap E,(ev) des composés MgsS et Mgo875X0125S (X=Cd, Zn) obtenues en utilisant 
l'approximation WC-GGA et comparées à d'autres valeurs théoriques et expérimentales disponibles 


dans la littérature. 
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Autres valeurs . . 
La direction 


théoriques de (a) 


du gap 


Exp{2] 5.19 Notre valeur 4.043, 
Notre valeur 79.21 
Exp[4] 5.201 


Expl4] Ref[10]4.01 
Ref[5]5.189, 


79.8 + 0.37, 76.0 + 
REf[4]5.18, Ref[8]4.248 


0.13 
Ref[6] 5.16, 
Ref[6], 81 
Ref[7]5.142 
Ref[9] 77.7 
Ref[8]5.234 


Figure IV. 1. Structures cristallines des composés M£sS (a), Mg 875Cdo 1258 (b) et Mg 87525 1258 (C) 


obtenues en utilisant l'approximation GGA-WC. 


IV. 3. Propriétés élastiques et mécaniques : 


Pour déterminer les caractéristiques mécaniques des composés Mg 375X91259 (X = Cd, Zn) 
dans la structure cubique RS, il suffit de calculer leurs constantes élastiques C;;, C;, et C44 via 
la méthode de Thomas Charpin, implémentée dans le code WIEN2K [11]. Ces paramètres 
fournissent des informations sur la manière dont un matériau retrouve sa forme initiale après 


la suppression d’une contrainte qui a subi. 


Les valeurs de C;,, C;, et C,, de nos composés sont présentées dans le tableau. IV. 2 et sont 
comparées à d'autres valeurs théoriques de la littérature. Les valeurs de C,, des composés sont 


significativement élevées, indiquant leur grande résistance aux forces de compression 
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unidirectionnelles, celles de C;, et C4, conduisent au module de cisaillement. Enfin, les 


valeurs des trois paramètres satisfont aux critères de stabilité de Born-Huang [12] suivantes : 


C0 C0 (6. - ce.) > 0,et (c., + 2c;.) > 0 


Tableau IV. 2. Valeurs des constantes élastiques des composés M£sS et Mg 875X0 1258 (X = Cd, Zn), 


calculées en utilisant l'approximation WC-GGA. 


M£g,.s75Cdo1258 79.75 35.39 38.26 
Mpgos752M5.125S 104.35 44.70 49.26 


Par conséquent, les composés Mg; 375X91259 (X = Cd, Zn) sont mécaniquement stables dans la 


structure cubique non magnétique. En utilisant les constantes élastiques, nous pouvons 
accéder à d'autres paramètres mécaniques tels que le module d'Young (E), le module de 
compressibilité (B), le coefficient de Poisson (v ), le module d'onde de cisaillement C,, le 
module de cisaillement Gx et le facteur d'anisotropie A via l'approximation de 


Voigt-Reuss-Hill [14,15,16]. Ces paramètres sont répertoriés dans le tableau IV. 3. 


La rigidité des matériaux peut être prédite connaissant les valeurs de leurs modules de Young 
E et de cisaillement G;, car on sait que lorsque E augmente, la rigidité du matériau augmente 
également. En fait, si l'on compare le MgsS aux composés Mg, 75Cdo 1258 et Mg 875200 1258, On 
peut facilement constater que le MgS (E= 196.38 GPa) est nettement plus rigide. Selon le 
rapport de Pugh B/Gx [17] et le coefficient de Poisson v, ces semi-conducteurs sont 


considérés comme durs. Il est connu que lorsque B/G;; est supérieur à 1.75, le matériau est dit 
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ductile ; sinon, il est considéré comme fragile. Nos composés Mg, 85X0 1258 (X = Cd, Zn) sont 


fragiles car les inégalités B/G; < 1.75 et v < 0.26 sont satisfaites. 


La stabilité structurale d'un composé peut être connue si son facteur d'anisotropie (A) est 
déterminé [18]. Les matériaux élastiques sont dits isotropes lorsque leur valeur de A est égale 
à 1 ; cependant, les écarts par rapport à l'unité de ce facteur se traduisent par une anisotropie 
élastique. Selon ce critère, le Mes et le Mgos75Zn5:S sont des matériaux isotropes, tandis 


que le Mg, 875Cdo 58 est anisotrope (A= 1.72). 


Nous avons calculé la température de Debye 8, et la température de fusion Le de nos 


composés, la connaissance de ces deux paramètres permet de remonter à leur chaleur 
spécifique, leur point de fusion, ainsi que leur expansion thermique. L'équation utilisée pour 


estimer 8, d’un composé est la suivante [19,20] : 


4, 

, sh 3n \ PV, LE 
D 1) Fe QE 
Fe 2 : + (IV. 2) 


1e __ [3B+4G 
mie v, =) 5 (IV. 3) 


Les symboles h, K3, N,, n, M et p représentent respectivement la constante de Planck, la 
constante de Boltzmann, le nombre d'Avogadro, le nombre d'atomes dans une seule cellule, le 


poids moléculaire et la densité du composé. V,,, v, v,, représentent respectivement la vitesse 


moyenne du son (m/s), la vitesse transversale du son (m/s) et la vitesse longitudinale du son 
(m/s). 


La température de fusion des composés de structure cubique est estimée à partir de la formule 


suivante [21] : 
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T. (K) = 553 + 5.91 (C11) + 300 (IV. 4) 


Le tableau IV. 3. Contient les valeurs calculées des grandeurs p, Vos vs, S) n'AT, des 


composés Mgos75X0125S (X = Cd, Zn). Le MgsS a la densité p (g/cm*)*10* la plus élevée, en 
second lieu vient le composé Mg: s75Cd5.5S et enfin le composé Mg 375Z/n0125S qui présente la 


densité la plus faible. Concernant les vitesses du son transversale v, (m/s) et longitudinale v, 


(m/s), qui représentent respectivement les vitesses des vibrations se déplaçant 
perpendiculairement ou dans la même direction de propagation dans le matériau, celles du 
composé Mgsg75/h55S sont les plus élevées et leurs valeurs respectives sont égales à 
9908.74 m/s et 16160.53 m/s. La vitesse moyenne du son calculée des composés 
Mgos875Cdo125S et Mgo875Z09.5S, ayant une liaison plus forte, montre des vitesses plus élevées 


que celle du Mgs pur. 


Rappelons qu’il est toujours possible de détecter un mode de vibration des phonons dans 


l'intervalle de températures délimité par les températures 8, etT, tel que 8, <T<T,, cet 


intervalle peut être exploré dans une étude vibrationnelle des composés. 


Tableau IV. 3. Valeurs calculées du facteur d'anisotropie élastique A, du module d'onde de 
cisaillement Cs, du module de cisaillement G;; (GPa), du module d'Young E (GPa), du module de 
compressibilité B (GPa), du rapport de Poisson © et du rapport B/G4, de la densité p du matériau, de la 
vitesse moyenne du son (m/s), de la vitesse transversale du son v, (m/s), de la vitesse longitudinale du 


son 7, (m/s), de la température de Debye 6, (K) et de la température de fusion T,, (K). 


__ 87.39 22.18 47.82 
C 


89.91 30.74 48.68 
Gy(GPa) 

196.31 76.59 120.51 
E (GPa) 

80.13 50.17 64.58 
B (GPa) 
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0.24 0.24 
V 
0.89 1.63 1.57 
B/Gy 


o (g/em*)x10° 179758.64 52537.92 49580.80 


7100 7649.19 9908.74 
vt (m/s) 


10548.29 13172.18 16160.53 
vl (m/s) 
7725.70 8487.69 10937.65 
V,, (m/s) 
772.90 530.96 691.87 
0» (K) 
1715.20 + 300 1024.32 + 300 1169.70+ 300 
Ty Œ) 


IV. 4. Propriétés électroniques : 


IV. 4.1. Bandes d'énergie : 


Le calcul des bandes d'énergie des composés MgS, Mgos75Cdo5S et MgosrsZNp15S à été 
exécuté en utilisant l'approximation TB-mBJ. Les bandes sont représentées le long des points 
de haute symétrie dans la première zone de Brillouin (figures IV. 2. 3. et 4). Le Mes pur 
présente un comportement semi-conducteur avec une large bande interdite indirecte d'une 
valeur égale à E,=4.043 ev dans la direction I°-X [22]. D'autre part, et sous l’effet du dopage, 
les composés respectifs Mg, s75Cdo5S et Mgos75Zo 58 sont des semi-conducteurs ayant des 
bandes interdites directes dans la direction '-[° de valeurs égales à 3.094 eV et 2.82 eV (voir 


tableau IV.1). 
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La bande interdite directe est une caractéristique recherchée dans les matériaux 


semi-conducteurs, car elle les rend très adaptables aux composants optoélectroniques tels que 


les lampes LED), les cellules solaires et les lasers [23,24]. 


MgsS-RS 


Energy(eV) 


e 


\ L r X  W Kk 


Figure IV. 2. Structure de bandes d’énergie du M£gs, calculée par l'approximation TB-mBJ. 
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Figure IV. 3. Structure de bandes d’énergie du Mg 375Cdo 1258, calculée par l'approximation TB-mBJ. 
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Figure IV. 4. Structure de bandes d’énergie du Mg 4752958, calculée par l'approximation TB-mBJ. 
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IV. 4. 2. La densité d'états : 


Les composés respectifs MgS, Mg s75Cd5 158 et Mg 875200 158, calculées via l'approximation 
TB-mBJ. La figure IV. 5 montre une densité électronique élevée vers le haut de la bande de 
valence près du niveau de Fermi, correspondant aux états électroniques de type p provenant 
du soufre (S) 

On note la présence des états s-Mg, p-Meg, s-S et p-S dans la bande de conduction. Les 
courbes TDOS des composés Mg 875Cdo125S et Mgos75Zo.15S présentées dans les figures IV. 
6. et 7 révèlent le caractère semi-conducteur à cause de l’absence d’états électroniques au 
niveau de Fermi. Les courbes TDOS des figures IV. 6. (a) et 7. (a) et les courbes PDOS des 
figures IV. 6. (b, c et d) et 7. (b, c et d) respectivement ont montré une contribution 
significative des états Cd-4d- t2g et Zn-3d-eg dans l'intervalle de -7 eV à -7.7 eV, et de -5 eV 
à -6 eV dans la bande de valence, ainsi que la présence des états S-p et une faible densité 
propres aux états Mg-p, Cd-4d-eg et Zn-3d-eg dans la zone située entre -4.2 eV et E;. Dans la 
bande de conduction, c’est les états Mg-s, Mg-p situés entre 3 eV et 16 eV qui sont 


dominants, en plus d'une faible présence d'états s-Cd et s-Zn entre 3 eV et 5 eV. 
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Figure IV. 5. Courbes de densités d'états totale et partielles du composé binaire M£gs, calculées en 


utilisant l'approximation TB-mBJ. 
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Figure IV. 6. Courbes de densité d'états totale (a) et partielles (b), (c), (d), du composé ternaire 


Mgos75Cdo5, calculées en utilisant l'approximation TB-mBJ. 
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Figure IV. 7. Courbes de densité d'états totale (a) et partielles (b), (c), (d), du composé ternaire 


Mgos75Z00 5, calculées en utilisant l'approximation TB-mBJ. 


IV. 5. Propriétés optiques : 


Les propriétés optiques d'un matériau donné constituent sa réponse à l’excitation 
électromagnétique, cette réponse contient des informations sur le matériau, principalement 
liées aux transitions électroniques qui peuvent avoir lieu quand il absorbe la lumière. 
Différents paramètres optiques peuvent être étudiés pour décrire optiquement un matériau et 
extraire ces caractéristiques, parmi lesquels nous citons la fonction diélectrique £(w), l'indice 
de réfraction n(w), le coefficient d'extinction k(w), la conductivité optique o(w) et le 
coefficient d'absorption a(w). Ces paramètres sont cruciaux à connaître pour accéder aux 


caractéristiques des matériaux. La fonction diélectrique £(w) est donnée par £(@) = &,(©) Hi 


e(&). 
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La partie réelle &,(w) est obtenue à partir de la partie imaginaire &,(w) par la transformation de 
Kramers-Kronig [7, 25]. Les quantités optiques &.(w), s(&), n(œ) et k(w) sont données par les 


expressions suivantes : 


e Q) = 1 + 2p [do (IV. 5) 
0 2 2. 
pi Po) à 
e (&@) =— +p J Fi dw (IV. 6) 
n(w) = eco) + Ru) (IV. 7) 
k(w) = 4 le(w)| — R)*) (IV. 8) 


IV. 5.1 Parties réelle £,(w) et imaginaire £,(©) de la fonction diélectrique : 


Nous avons obtenu à partir des bandes d'énergie les valeurs du gap énergétique (Eg) des 
composés MgeS, Mgos75Cdo5S et Mgog75Z55S, qui sont respectivement de 4,043 eV, 3,094 
eV et 2,82 eV. Nous avons relevé les valeurs de &, (0), qui désignent la partie réelle de la 


fonction diélectrique statique (valeur à énergie nulle), des courbes présentées dans la figure 
IV.8, elles sont égales à 4.64, 4.87 et 4.94 pour les composés respectifs MgS, Mg 875Cdo 1258 


et Mgos5Zno1sS. On remarque que les valeurs de E, diminuent au fur et à mesure que la 
valeur de la constante diélectrique £, (0) augmente, conformément au modèle de Penn [26] 


selon lequel E, et £,(0) sont liés par l'équation (9) ci-dessous : 


2 
€ (0)= 1 + Fe) (IV. 9) 


g 


Les principaux pics de E (w) dans la figure. IV. 8 sont principalement situés autour de 6.35 


eV, 5.76 eV et 6.20 eV, ils proviennent de transitions électroniques impliquant les états 
électroniques S-3p et Mg-s et Mg-p. Dans les trois composés, l'amplitude de ces pics 
augmente avec l'énergie de l'onde et atteint un maximum entre 9 eV et 10 eV. Au-delà de 10 


eV, € ,(@) devient presque nul. 
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Figure IV. 8. Parties réelles E,(@) de la fonction diélectrique des composés MgsS et Mg 375 X 91259 (X 
= Cd, Zn), calculées par l'approximation TB-mBJ. 
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Figure IV. 9. Parties imaginaires £, (w) de la fonction diélectrique des composés M£s et 


Mgo875X01259 (X = Cd, Zn), calculées par l'approximation TB-mBJ. 


La figure IV. 9. exprime l’évolution de la fonction € .(@) en fonction de l'énergie d'onde dans 


l'intervalle entre O0 et 40 eV. Nous avons noté le pic le plus significatif dans le spectre 
ultraviolet aux énergies 7.03, 6.99, 7.08 eV pour les composés respectifs MgS, Mg, 875Cdo 1258 


et Mgog75Zho:S. Ces observations notées sur les deux spectres E,(w) des composés 
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indiquent des transitions électroniques significatives dans cette région. L'analyse des courbes 
TDOS et PDOS (les figures IV. 5. 6. et 7.) des composés permet de comprendre que le pic 
principal provient des transitions entre les états électroniques de type S-3p situés en haut de la 
bande de valence et les états électroniques inoccupés de type Mg-3s et Mg-2p situés en bas de 
la bande de conduction. Pour les composés ternaires Mg, 375Cd5 158 et Mgog75Z9 1259, nous 
avons observé la présence d'états électroniques dans la gamme de 2 eV à 6 eV. 

Ces pics correspondent à des transitions inter-bandes entre des états électroniques occupés de 
types S-3p, Cd-4d-eg / Zn-3d-eg situés près du niveau de Fermi dans la bande de valence et 
des états électroniques inoccupés de types Cd- Ss / Zn-4s situés dans la bande de conduction. 
Entre 10 eV et 40 eV, les courbes des composés MgS, Mg 375Cdo 158 et Mgo 87520 58 ont des 


allures très similaires. 


IV. 5.2. Indice de réfraction n (o) et indice d'extinction Kk (o) : 


L'indice de réfraction n (œ) et l'indice d'extinction k (æ) sont essentiels pour comprendre 
comment la lumière se propage dans un matériau. La figure IV. 10. (a) montre la variation de 
n(æ) en fonction de l'énergie entre O et 40 eV des composés MeS, Mgoe5Cdo5S et 
Mgos75Zno 8. Les courbes ont des comportements identiques et présentent un pic intense 
situé dans la région de la lumière ultraviolette entre 4 eV et 7.4 eV. Les positions de ce pic 
dans les composés respectifs MgS, Mgo375Cdo125S et Mgo875Zn0 1258 Sont les suivantes : 6.45 
eV, 5.85 eV et 6.27 eV. Au-delà de 7 eV, n(œ) diminue, sa valeur la plus faible dans les trois 
composés est n(o) = 0.25, puis son amplitude continue à augmenter légèrement. Les valeurs 
statistiques n (0) des composés MgS, Mgos75Cdo5S et Mgog75Zh0125S Sont respectivement 
égales à 2.15, 2.20, 2.22. 

La variation de l'indice d'extinction K(®) en fonction de l'énergie du M£gS pur et des 
composés ternaires sont montrées dans la figure IV. 10. (b). Cette quantité, dont la forme 


ressemble à celle de £ .(@), varie significativement avec l'énergie absorbée de l'onde. Le pic 


intense de K(o) est situé à 7.3 eV, 7.11 eV, et 7.30 eV dans les composés respectifs Mgs, 
Mgos75Cdo5S, et Mg 875200 8. De plus, les courbes K(w) de tous les composés diminuent 


significativement aux énergies supérieures à 12 eV. 
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Figure IV. 10. (a) Indice de réfraction n (&) et (b) Indice d'extinction K (@), des composés Mes et 
Mgos875X01259 (X = Cd, Zn), calculés par l'approximation TB-mBJ. 


IV. 5.3. Conductivité optique et coefficient d'absorption : 


D'après la figure IV 11. (a), la courbe de conductivité optique © (w) du MgsS débute à 5.4 eV 
et celles des composés ternaires Mg, 475Cdo 1258 et Mg 875205 1258 à 4.8 eV. Le pic intense est 
observé aux énergies approximatives égales à 7.00 eV, 7.11 eV et 6.99 eV dans les composés 
respectifs Mes, le Mg, s75Cdo 58 et le Mg 875200 58, celui du M£gs est le plus prédominant et 
sa valeur est de 15331.5 (Ohm. cm) '. Aux énergies supérieures à 9.23 eV, l'amplitude de o(o) 
diminue dans les trois composés. 

Le coefficient d'absorption a(w) fournit des informations sur la capacité d'un matériau à 
absorber la lumière. La figure IV. 11. (b) montre la variation de a(w) des composés en 
fonction de l’énergie de l’onde dans l’intervalle entre 0 et 40 eV. Les calculs de a(w) sont 


effectués via l'approximation TB-mBJ. Le principal pic de a(w) apparait dans la région de la 
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lumière UV entre 6 eV et 12 eV avec une amplitude égale à 227.63 (10*.cm-1) pour Mgs, 
208.55 10*.cm' pour Mg, g75Cdo 58 et 215.13 10*.cm'! pour Mg, 7529 158. Au-delà de 16 eV, 


l'amplitude de a(w) diminue, indiquant une faible absorption de la lumière par les composés. 
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Figure IV. 11. (a) Conductivité optique o(@) et (b) coefficient d'absorption a(®) des composés M£s et 
Mgos875X0.1259 (X = Cd, Zn), calculés en utilisant l'approximation TB-mBJ. 


IV. 6. Propriétés thermoélectriques : 


Les propriétés thermoélectriques déterminent collectivement l'efficacité d'un matériau à 
exploiter la chaleur perdue ou les gradients de température pour des applications pratiques 
telles que la production d'énergie, le refroidissement des ordinateurs et la réfrigération [27]. 
Les paramètres thermoélectriques à investiguer pour se prononcer sur les performances des 
matériaux à très haute température sont les suivants : un coefficient de Seebeck élevé, une 
excellente conductivité électrique et une conductivité thermique minimale. 

Pour les semiconducteurs, ces propriétés sont excellentes, et c’est pour cette raison que nous 
avons vu important d’étudier les propriétés thermoélectriques de nos composés. Dans cette 


section, nous calculons les paramètres thermoélectriques des composés Mg35X0 1258 (X = 
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Cd, Zn) à savoir, le Coefficient de Seebeck (S) (figure IV. 12), la Conductivité électrique (6) 
(figure IV.13), la Conductivité thermique électronique (K) (figure IV.14), la Conductivité 


thermique du réseau (k, ) (figure IV.15), la figure de mérite ZT (figure IV.16) et le facteur de 
puissance (PF) (figure IV.17). Ces paramètres sont investigués à l'aide du package Boltztrap 
[28] basé sur un algorithme semi-classique avancé du transport de Boltzmann pour étudier les 
propriétés thermoélectriques en exploitant la densité locale des états et la structure des bandes. 
Nous calculons le facteur de mérite ZT selon la relation suivante [29] : 


2 
ZT = Pia (IV.10) 

Où PF = S?c est le facteur de puissance [30], et k, + k, sont la somme de la conductivité 
thermique électronique k, et de la Conductivité thermique du réseau k, . La conductivité 
thermique k, est calculée par le package Boltztrap comme k, = LTo, qui est basé sur la loi de 
Wiedemann-Franz [31], où L est le nombre de Lorentz. Cependant, la Conductivité thermique 
du réseau k, n'est pas calculée, mais son calcul est nécessaire et s’effectue à partir de la 


relation de Slack suivante [32] : 


k = 4. (IV. 11) 
l Tr 


Où M , 0» V,T, n, y, représentent respectivement la masse atomique moyenne, la température 
de Debye, le volume de la maille, la différence de valeur de température, le nombre d'atomes 
dans la maille élémentaire, et y est le paramètre de Grüneisen calculé par l'expression 


proposée par Julian [33] : 


ÿ= — (IV. 12) 


où v, v, sont respectivement les vitesses longitudinale et transversale du son. 


La constante À dépend du paramètre de Grüneisen, soit : 


2.43x107 
À — ER (IV. 13) 


F: 
Y y 


Le temps de relaxation utilisé dans le code Boltztrap proposé par Ashcroft et Mermin [341] est 


constant et sa valeur est t = 0,8 x 10°*s. 
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Le coefficient Seebeck (S) est une mesure de la capacité d'un matériau à générer un potentiel 
électrique en réponse à un gradient de température. La figure IV. 12. Présente les courbes de 
(S) des composés M£gS, Mgog75X01258 (X = Cd, Zn) dans la gamme de température variant 
entre 50 et 1200 K. La courbe (S) du composé binaire MgS diminue avec l'augmentation de 
la température jusqu’à 450 K, elle varie ensuite d’une manière monotone jusqu’à 1200 K et 
atteint une valeur minimale de 240.7.10° [uV/K]. Concernant le coefficient Seebeck des 
composés Mg 875Cdo125S et Mgog7s Zho5 (S), il ne possède pas la même allure que celle de 
Mes pur, on note une amplitude plus importante du composé dopé au cadmium (Cd) 
comparée à celle du composé dopé au zinc (Zn). Les deux courbes (S) des composés ternaires 
Mg-Cd-S et Mg-Zn-S présentent respectivement des maximums de valeurs égales à 
269.41.10% [uV/K ] à T-201.6 K et à 209.41.10° [uV/K ] à T-250.96 K. Au-delà de ces deux 
maximums, les courbes (S) des deux composés ternaires suivent une allure décroissante et 
atteignent à T—1200 K les valeurs suivantes : 202.71.10° [u V/K] pour le composé dopé au Cd 
et 195.53.10° [u V/K] pour le composé dopé au Zn. Rappelons que les valeurs positives de (S) 


indiquent que les composés ternaires sont des semi-conducteurs de type p. 
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Figure IV. 12. Coefficient de Seebeck (S) des composés Mes et Mg, 375X 91258 (X = Cd, Zn) entre 50 K 


et 1200 K. 


Les courbes de conductivité électrique (6) des composés Mg et Mg, 375X 01258 (X = Cd, Zn) 


sont présentées dans la figure IV. 13, entre 50 K et 1200 K, elles suivent une allure croissante 
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avec l’augmentation de la température. Cependant, il est à noter que la conductivité électrique 
6 du composé dopé au Zn est plus importante que celles des deux autres composés, elle atteint 
la valeur 2.5 .10° (Q. m) * à T= 1100 K. Rappelons qu’une conductivité électrique élevée des 
composés Mgs75X01259 (X — Cd, Zn) est essentielle pour une conversion efficace de la 
chaleur en électricité, elle permet la circulation des électrons, contribuant ainsi à la génération 


d'une tension thermoélectrique. 
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Figure IV. 13. Conductivité électrique (6) des composés M£sS et Mg 875 X0 1259 (X = Cd, Zn) entre 50 
Kelvin et 1200 Kelvin. 


Les figures IV. 14, 15, représentent respectivement la conductivité thermique du réseau k, et la 
conductivité thermique électronique k,, k, représente la capacité du réseau cristallin à conduire 
la chaleur associée aux vibrations des phonons. La figure IV. 14. illustre la variation de k, en 
fonction de la température des composés Mges et Mg, 575 X0 1258 (X = Cd, Zn). Comme il est vu 
sur la figure IV. 14. k, est une fonction constamment décroissance avec l'élévation de la 
température et ceci pour tous les composés, les courbes se stabilisent aux valeurs de k, 


suivantes : 1.51.10° [W/m.K] pour le MgS , 445.10 [W/m.K] pour Mgss5Cdo5S, et 
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8.65.10TW/m.K] pour le Mgog75Zo 58. Rappelons qu’il est souhaitable de minimiser la 


conduction de la chaleur à travers le composé pour un meilleur rendement thermoélectrique. 
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Figure IV. 14. Conductivité thermique du réseau (K;) des composés Mes et Mo, 875X 01259 (X = Cd, Zn) 


entre 50 K et 1200 K. 


La conductivité thermique électronique k, provient de la dynamique des porteurs de charge et 
de la capacité des électrons à transférer la chaleur sous l'influence de la température. Dans les 
métaux et les semi-conducteurs, les électrons contribuent de manière significative à la 
conduction thermique. La figure IV. 15 Présente la variation de k, en fonction de la 
température des composés M£gS et Mg535X01258 (X = Cd, Zn), elle montre une croissance 
continue de k, quand la température augmente entre 50 K et 1200 K. Pour le Mgs pur, les 
valeurs de k, varient entre 6.44.10"? [W/ m.K.s] à 100 K et 2.25.10% [W/ m.K.s] à 1200 K, 
concernant les composés ternaires, la valeur de k, varie de 1.44.10*[W/m.K.s] à 100 K 
jusqu'à 1.62.10% [W/ m.K.s] à 1200 K pour Mgsg5Cdos S et de 3.06.10'* [W/ m.K.s] à 100 


K jusqu'à 2.19.10 *[W/ m.K.s] à 1200 K pour Mgo 875 Zno 58. 
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Figure IV. 15. Conductivité thermique électronique (Kk.) des composés Mgs et Mg, 85X0 1258 (X = Cd, 


Zn) entre 50 K et 1200 K. 


La figure IV. 16 exprime la variation du facteur de puissance (PF) des composés M£gs et 


Mgos875X01258 (X = Cd, Zn) dans l’intervalle de températures allant de 50 à 1200 K. 


Le facteur de puissance dépendant de la température est essentiel pour adapter les matériaux à 


des performances thermoélectriques optimales. Il s'agit d'un compromis délicat qui vise à 


maximiser le coefficient Seebeck et la conductivité électrique et en même temps vise à 


minimiser la conductivité thermique et tout ça dans le but de récolter et de convertir 


efficacement l'énergie. On observe clairement sur la figure IV. 16, pour tous les composés, 


que l’amplitude de PF augmente avec la température sur tout l’intervalle de température 


exploré ; néanmoins, elle est plus importante entre 250 K et 650 K pour les composés 


ternaires. 
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Figure IV. 16. Facteur de puissance des composés MgsS et Mg, 875X01259 (X = Cd, Zn) entre 50 K et 


1200 K. 


Le facteur de mérite ZT représente un autre paramètre essentiel utilisé pour évaluer les 
performances thermoélectriques des matériaux. Les courbes ZT en fonction de la température 
des composés sont présentées sur la figure IV. 17, entre 50 K et 1200 K, notons d’abord qu’il 
y a une ressemblance entre les courbes ZT et les courbes du coefficient de Seebeck (S). Pour 
le binaire MgS, la courbe ZT subit une décroissance monotone avec l’augmentation de la 
température et arrive à la fin du spectre à la valeur ZT=— 6.01.10°"; inversement, la valeur de 
ZT des deux composés dopés augmente avec l'augmentation de la température depuis le début 
du spectre et passe par un maximum de 5.93.10°% à 451.09 K pour Mgos75Zo 8 et par un 
maximum de 6.47.10°% à 351.02 K pour Mgog75Cd5 58. Notons bien que la valeur ZT du 


composé dopé au cadmium (Cd) est meilleure que celle du composé dopé au zinc (Zn). 
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Figure IV 17. Figure de mérite (ZT) des composés M£sS et Mg, 875X0 1258 (X = Cd, Zn) entre 50 K et 
1200 K. 
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Chapitre V 


Etude des propriétés électroniques et optiques du 
semi-conducteur M£,375X 91259 (X =Co, Ni) dans la phase 
antiferromagnétique. 
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Introduction : 


Dans cette section, l'investigation se focalise sur les composés M£sS dopés avec les métaux de 
transition cobalt (Co) et nickel (Ni). L'analyse est réalisée sur les composés Mg, 85X0 1258 (X 


= Co, Ni) dans la structure antiferromagnétique rock salt (RS). 


Les propriétés des composés étudiés sont de types structuraux, électroniques et optiques. 


Dans ce type de composés, les moments magnétiques des porteurs adjacents s'alignent dans 
des directions opposées, entraînant une annulation de l'aimantation macroscopique. Ces 
propriétés rendent ces matériaux particulièrement intéressants pour leur potentiel dans le 
stockage de données et la magnéto-électronique, et ils sont utilisés dans la conception de 
capteurs magnétiques très sensibles pour des applications telles que la mesure de champs 


magnétiques faibles ou la détection de mouvements [1,2]. 


V. 1. Détails des calculs : 


Les caractéristiques structurales, électroniques et optiques des composés Mg 375X 91259 (X = 
Co, Ni) dans la structure rock salt antiferromagnétique sont déterminées à l'aide de la méthode 
des ondes planes augmentées linéarisées à plein potentiel (FP-LAPW), mise en œuvre dans le 
logiciel de calcul Wien2K et basée sur la théorie de la fonctionnelle de la densité (DFT). 

Les propriétés structurales sont déterminées à l'aide de l'approximation du gradient généralisé 
proposée par Wu et Cohen (WC-GGA), tandis que les propriétés électroniques et optiques 
sont calculées par l'approximation TB-mBJ. Pour effectuer les calculs des propriétés des 
composés ternaires Mgog75X012589 (X = Co, Ni) dans la phase antiferromagnétique nous 
remplaçons deux atomes de Mg par les éléments de transition Co et Ni tout en respectant un 
alignement anti-parallèl des moments magnétiques des porteurs. La maille de réseau cubique 


de chaque composé est créée sous forme d'une supercellule bec a 2X2X2 contenants 32 
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atomes, où les atomes X = Co, Ni occupent dans la maille cubique les sommets (0,0,0), et les 
milieux des faces (0.5, 0, 0) (Voir la figure V.1). 

La valeur choisie pour le rayon muffin tin de la sphère atomique est telle que Rur.K,,, = 8, où 
Kax est la valeur maximale de l'onde dans la zone de Brillouin. Les valeurs du rayon muffin 
tin (Rur, des atomes Co/Ni, Mg et S sont respectivement égales 2.5, 2.3, 2.04 (a.u). La valeur 
maximale de l'onde partielle à l'intérieur de la sphère atomique est 1, = 10. De plus, la 
densité de charge de Fourier est étendue jusqu'à G... = 16 pour capturer avec précision les 
caractéristiques essentielles, en utilisant 1000 points -k pour une analyse complète des 
propriétés structurales et électroniques et 3000 points-k pour le calcul des propriétés optiques. 
Les critères de convergences de l'énergie totale et de la densité de charge des électrons sont 


respectivement fixés à 10 *Ry et à 10° Ry. 


V. 2. Propriétés structurales : 


Le paramètre de réseau a (À) et le module de compressibilité B (GPa) ainsi que sa dérivée B' 
sont calculés en utilisant les paramètres ajustés à partir de l'équation d'état de 
Birch-Murnaghan [4]. 

Les valeurs de a (À) et de B (GPa) du MgsS pur sont en parfait accord avec les valeurs 
expérimentales et théoriques de la littérature [5,10], elles sont présentées dans le tableau V. 1. 
Concernant les composés Mgos75Co01sS et MéossNiowsS, les valeurs de a (À) sont 
respectivement égales à 5.16 (À) et 5.15 (À) et sont relativement faibles par rapport à celle de 
M£gsS pur et ceci en raison de la différence entre les rayons ioniques de ce dernier et des 


atomes Co, Ni. 


La stabilité magnétique des composés Mg, 37:X9 1258 (X = Co, Ni), est obtenue en optimisant 


l'énergie des composés en fonction du volume de la maille élémentaire, les calculs ont été 
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faits dans la phase ferromagnétique (FM) et dans la phase antiferromagnétique (AFM), le 
principe consiste à chercher l’énergie la plus faible des deux phases pour décider laquelle des 


deux est la plus stable. Les résultats sont rangés dans le tableau V.1 où 


ÂEË = AEyry — AEry est la différence entre les énergies des phases antiferromagnétique et 
ferromagnétique. 

Si AE < 0, cela signifie que la phase AFM est la plus stable pour les composés ; si AE > 0, 
c’est la phase FM qui est la plus stable. Les valeurs du tableau V.1 montrent des valeurs 


négatives pour ÂE indiquant, ainsi, la stabilité des deux composés dans la phase AFM. 


Tableau V. 1 Valeurs calculées du paramètre de réseau a (À), du module de compressibilité B (GPa) 
et de sa dérivée première B', des énergies minimales des phases FM (E;u(eV)) et AFM (Exrm(eV)) et 
de la différence AE entre les énergies des phases des composés Mgss75X91258 (X = Co, Ni). Les 


calculs sont effectués en utilisant l'approximation WC-GGA et sont comparées aux valeurs théoriques 


et expérimentales publiées antérieurement. 


Composé Autres Autres 
a(À) calculs calculs de 
de a(À) 


| 


Mgos875C00:125S 
—0.00 


Mgos7sNio.125S ; | 0001 


131 


Figure V. 1. Structure cristalline et arrangement des moments de spins dans la phase 


antiferromagnétique : (a) Mgos75C06125S et (b) Mg 875Nio 1258 : 


V. 3. Propriétés électroniques : 


V. 3.1. Structures de bandes : 


Les bandes d'énergie des composés Mg 75X 91258 (X = Co, Ni) ont été calculées le long des 
directions de haute symétrie (R-[°) dans la zone de Brillouin en utilisant la méthode TB-mBJ. 
Les figures V2 et V3 contiennent les bandes d’énergie des composés Mgo375C001258 et 
Mgos875Ni0125S et montrent leur caractère semi-conducteur dans le canal des spins majoritaires 


(f) et dans le canal des spins minoritaires (|). 


Le composé Mg,375C0015S présente un gap indirect dans la direction M-X dans les deux 
canaux de spins dont les valeurs sont égales à 1.15 eV dans le canal des spins majoritaires et 
1.10 eV dans le canal des spins minoritaires, la différence entre les deux valeurs et faible de 
l’ordre de 5%. Le composé Mg, 45Ni5.5S à aussi un gap indirect dans la direction M-T pour 
les deux populations de spins (up+down) et dont les valeurs sont égales à 1.46 eV et 1.48 eV 


dans le canal des spins majoritaires et dans le canal des spins minoritaires respectivement, 
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l’écart entre les deux valeurs n’est que de 2%. Les valeurs du gap énergétique des différents 


composés (binaire et ternaires) sont contenues dans le tableau V.2. 


En plus des valeurs du gap énergétique des composés, le tableau V2 contient également les 
valeurs de leurs moments magnétiques totaux MM' et locaux MM* relatifs aux éléments qui 
constituent les composés. En raison d’une distribution antiferromagnétique des moments 
(alignement anti-parallel), les moments totaux sont donc nuls. 


Tableau V. 2. Valeurs du gap énergétique (E,) et des moments magnétiques totaux et locaux des 


composés Mg, 875X 0.159 (X = Co, Ni), calculées en utilisant la méthode TB-mBJ. 


Composé | MM° | MM° | MM | MM |E,(eV)| EE," (eV) | E,“ (eV) 
(us) | (us) 5) (5) 


M£S 4.043, 
4.010! 
Mara CoueS Fo nu og HAE EE 5 
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Figure V. 2. Bandes polaris 


FigureV. 3. Structures de bandes polarisées en spin du composé Mg s5Nio 5, calculées en utilisant 
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V. 3. 2. Densité d’états électroniques : 


Les densités d'états polarisées en spins totales (TDOS) et partielles (PDOS) des composés 
Mgo875X0125S (X = Co, Ni) ont été calculées à l'aide de l'approximation TB-mBJ dans la 
structure rocksalt antiferromagnétique. Elles sont présentées dans les figures V.4 et V.S sur 
une plage d'énergie allant de -8 eV à 13 eV, pour les deux canaux de spins, majoritaire (up) et 


minoritaire (dn). 


Les courbes TDOS V4 (a), et V.S (a) des deux composés ternaires indiquent un caractère 
semi-conducteur, elles sont symétriques dans les deux canaux de spins du fait de 
l’anti-alignement des moments magnétiques des états 3d-Co et 3d-Ni causé par l’ordre 


antiferromagnétique. 


Les figure V.4 (b) et V.4 (c), V.S (b)et V.5 (c) , montrent les bandes de valence des composés 
ternaires dans l’intervalle de -8 eV à 13 eV. Dans cette région, le composé Mg, 47:5C05 1258 fait 
apparaitre les états Co-3d-t2g entre -4.87 et -0.3 eV et le composé Mg ,875Nio125S montre les 
états Ni-3d-t2g entre -4.4 et -0.5 eV. Sur ces deux intervalles existent aussi les états 


électroniques de types « p » et «s » appartenant aux éléments Mg ets. 


Les sommets des bandes de conduction des deux composés sont principalement occupés par 
les états électroniques de type « s » et & p », tandis que leurs bas contiennent les états d-t2g 
propres aux éléments Co et Ni, ainsi que les états p-Mg, localisés à proximité du niveau de 


Fermi. 
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Figure V. 4 Courbes de densités d'états totale et partielles du composé Mg; 75C00 1258, calculées par 


l'approximation TB-mBJ . 
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Figure. V. 5. Courbes de densités d'états totale et partielles du composé Mg 75Nio 58, calculées par 


l'approximation TB-mBJ. 
V. 4. Propriétés optiques : 


Les propriétés optiques d'un matériau décrivent sa réaction lorsqu'il est excité par une onde 
électromagnétique. Dans cette section, nous exposons les propriétés optiques des composés 
Mgos875X01258 (X= Co, Ni) en calculant les grandeurs optiques suivantes : la partie réelle £, (&) 
et la partie imaginaire £, (@) de la fonction diélectrique, l'indice de réfraction n (), l'indice 
d'extinction k (w), la conductivité optique o(o), et le coefficient d'absorption à (æ). Tous ces 
paramètres sont obtenus en utilisant les relations de Kramers-Kronig [13]. Les calculs sont 


réalisés à l'aide de l'approximation TB-mBJ. 
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V. 4. 1. Partie réelle £1(w) et partie imaginaire £2(w) de la fonction 


diélectrique : 


Les figures. V.6 et V.7 présentent respectivement la partie réelle £&,(w) et la partie imaginaire 
£æ) de la fonction diélectrique du composé binaire MgS et des composés ternaires 


Mgos875X0125S (X= Ni, Co) en fonction de l’énergie du photon w dans l’intervalle de 0 à 35 eV. 
La figure. V.6 montre la présence d’un pic intense centré aux environs de 6 eV dans les 
composés ternaires et à une énergie légèrement supérieure à 6 eV dans le composé binaire 
Mes. Un deuxième pic est observé aux basses énergies dans les composés ternaires, il est 
situé aux énergies 1.22 eV et 1.99 eV dans Mgo375C00 1258 et Mgog75Ni9 1258 respectivement et 
provient des transitions entre des états occupés dans la bande de valence des composés vers 
des états inoccupés dans leur bande de conduction dans les canaux de spins majoritaires (up) 


et minoritaires (down). 


Tableau V. 3. Amplitudes et positions des pics relevées des courbes £1 des composés Mgo 875X 0 1258 


(X = Co, Ni). 


Position du pic (eV) Amplitude du pic 


. 


Les courbes £,(w) des composés ternaires sont identiques (superposées) au-delà de 2 eV, la 


différence, bien que faible, est observée aux basses énergies dans l’intervalle 0 à 2 eV. 
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Au-delà de 9.5 eV, les trois courbes sont identiques en termes d’amplitude. Entre 9.5 eV et 19 
eV, le signe de £&,(w) est négatif témoignant de la métallicité des composés et de leurs pouvoirs 
à réfléchir la lumière incidente. Aux énergies E>19 eV, &,(w) devient positif et suit une 
croissance monotone jusqu’à la fin du spectre optique où la valeur de son amplitude 


approxime 1. 


Les valeurs statiques de la partie réelle de la fonction diélectrique s.(0) sont respectivement 


égales à 4.64, 5.71 et 5.46 pour les composés MgS, Mg, 375C00 1258 et Mgog75Ni91259. 
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Figure V. 6. Parties réelles £1 de la constante diélectrique des composés Mg et Mg685X0 1259 (X = Co, 


Ni), calculées par l'approximation TB-mBJ. 


La figure V. 7 présente la variation de la partie imaginaire (£,) de la fonction diélectrique pour 
les composés MgesS et Mgs5X0 1258 (X= Ni, Co). On observe un pic additionnel aux faibles 
énergies dans l'intervalle de la lumière visible entre 1.5 à 3 eV dans les composés 
Mgos875C00125S et Mg9875Nio.15S, il est dû aux transitions électroniques mettant en jeu les états 
Co-3d-t,, et Ni-3d-t,, dans le canal des spins minoritaires (dn) dans la bande de valence vers 
les états p-Mg vides dans la bande de conduction. Dans la région ultraviolette de 6 à 9 eV, on 


observe deux pics dont l’intensité est plus importante que celle du pic existant dans le visible. 
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Les amplitudes ainsi que les positions énergétiques des différents pics de £, des composés sont 


regroupées dans le tableau V. 4. 


Tableau V. 4. Amplitudes et positions des pics relevées des courbes £, des composés Mg, 55X 01259 


(X= Co, Ni). 


Position du pic (eV) Amplitude du pic 


Selon les courbes TDOS et PDOS (figures V. 4 et V. 5.), les pics observés sont attribués aux 
transitions entre les états électroniques de type S-3p, localisés dans la région supérieure de la 
bande de valence, et les états électroniques non occupés de type Mg-3s et Mg-2p, localisés 


dans la portion basse de la bande de conduction. On remarque que la variation de E,(@) de 


ces deux composés est similaire avec une très légère différence dans la région de l’infrarouge. 


À haute énergie (E > 9 eV), les courbes de £ .(@) diminuent et leurs amplitudes deviennent 


presque nulles au-delà de 26 eV. 
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Figure V. 7. Parties imaginaires €, de la constante diélectrique des composés M£sS et Mg 875 X 91258 (X 


= Co, Ni), calculées par l'approximation TB-mBJ. 


V. 4. 2. Indice de réfraction n (o) et indice d'extinction K (o) : 


L'indice de réfraction n (@) et l'indice d'extinction k (w), détaillant la transmission de la 
lumière à travers un milieu, sont illustrés dans la figure V. 8. (a). La variation de n (&œ) des 
composés, sur une grande partie de l’intervalle d’énergie exploré, ne présente pas de 
différences signifiantes, à l’exception de la région de la lumière infrarouge-visible qui montre 
un écart entre les formes et les amplitudes des courbes des composés où, en effet, on note la 
présence d’un pic dans les composés ternaires et son absence dans le composé binaire. Les 
positions ainsi que les amplitudes des pics de n (©) des composés sont rapportées dans le 


tableau V. 5. 
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Tableau V. 5. Amplitudes et positions des pics de n (œ) des composés Mg et Mg, 875X 91259 (X = Co, 


Ni). 


Composé Position du pic (eV) Amplitude du pic 
 S 


PE ME ME 


(l 


Le M£gsS possède à la valeur maximale de l'indice de réfraction à environ 6.4 eV qui vaut 3.7, 


dans les composés ternaires cette valeur est située approximativement à 6 eV et vaut 3.4. 


Le troisième pic observé vers 8 eV dans les trois composés à une intensité un peu plus faible 
que celle du pic qui le précède et vaut à peu près 2.25. Au-delà de 8 eV, n (œ) décroit en 
fonction de l’énergie jusqu’à 16 eV, ensuite il reprend une croissance lente jusqu’à la fin du 
spectre où son amplitude est approximativement égale à 0.8. Au zéro des énergies, les valeurs 
de l’indice de réfraction statique n (0) des composés MgS, Mg, 875C00 158 et Mgo 875 Ni 1259 
sont respectivement égales à 2.15, 2.39, 2.33, la valeur la plus élevée est celle du composé 
dopé au nickel (Ni), il est clair que le dopage du composé binaire par les métaux de transition 


a conduit à une augmentation de sa valeur de n (0). 


Les courbes de l'indice d'extinction k (®) des composés M£sS et Mg, 55X0 1258 (X = Co, Ni) 
sont illustrés dans la figure. V. 8. (b), elles présentent deux pics intenses dans la région 


ultraviolette aux positions énergétiques contenues dans l’intervalle de 7 eV et 9.2 eV, leurs 
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amplitudes sont égales à 2.5 et 2.6 pour Mes, 2.25 et 2.29 pour Mg, 475C05158 et 2.32 et 2.22 
pour Mgs375Ni0 159. Aux basses énergies, dans la région du visible, on note l’existence d’un 


pic dans les composés ternaires d'amplitudes approximant la valeur 0.5. 
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Figure V. 8. (a) Indice de réfraction n (&) et (b) Coefficient d'extinction K (œ), des composés M£s et 
Mgos875X01259 (X = Co, Ni), calculés par l'approximation TB-mBJ. 


V. 4.3. Conductivité optique et coefficient d'absorption : 


Les courbes de conductivité optique o(w) des composés MgsS et Mg, 75X0 1258 (X = Co, Ni) 
sont illustrées dans la figure V. 9. (a), elles sont tracées dans l’intervalle entre 0 et 35 eV et 
ont un profil similaire à celui de K (œ). Une forte conductivité optique est enregistrée dans la 
région de l’ultraviolet avec l’apparition de deux pics intenses aux positions situées 
approximativement à 7 eV et 9 eV. Les courbes des composés ternaires ont des profils 
similaires et sont presque superposées avec une légère différence insignifiante observée entre 


3 eV et 4 eV où l’amplitude du composé contenant le nickel est un peu plus grande que celle 
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du composé contenant le cobalt. Le deuxième pic relevé aux faibles énergies entre 2 eV et 3 
eV est spécifique aux composés ternaires et met en jeu les transitions électroniques 
impliquant les états 3d de l’élément de transition avec lequel le M£sS a été dopé, ce pic est 
inexistant dans la courbe o(œ) du composé binaire. Les positions et les intensités des pics de 


la fonction o(&) des trois composés sont montrées dans le tableau V. 6. 


Tableau V. 6. Amplitudes et positions des pics de o (w) des composés MgsS et Mgos75X01258 (X = Co, 


Ni). 


Composé Position du pic (eV) 


Au-delà de 9 eV, les trois courbes o(@) sont pratiquement supperposées et leurs amplitudes 
diminuent avec l’accroissement de l’énergie jusqu’aux grandes énergies. La décroissance de 
6(w) des trois composés est lente à partir de 16 eV alors qu’elle est rapide dans l’intervalle de 
9 eV à 16 eV. La valeur de c(w) est faible aux grandes énergies, elle est égale à 252.63 


[Q.cm] * à la fin du spectre dans les trois composés. 


La figure V. 9 illustre la variation du coefficient d'absorption a(w) des trois composés dans 
l'intervalle de 0 à 35 eV. Une forte absorption est enregistrée dans la région de l’ultraviolet, 
l’intensité maximale est atteinte aux environs de 9 eV. Comme pour c(w), les courbes a(@) ne 
présentent pas de différences signifiantes quand on passe du composé binaire vers les 
composés ternaires ni en termes de profil ni en termes d'amplitude excepté dans la région du 


visible où apparait un pic supplémentaire de faible amplitude et dont l’origine est liée à 
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l’élément de transition (Co, Ni) étant donné qu’il est absent sur la courbe de la matrice hôte 
semi-conductrice (M£gS). Les positions des pics ainsi que leurs intensités sont affichées dans 


le tableau V. 7. 


Tableau V. 7. Amplitudes et positions des pics de a(w) des composés MgsS et Mg, 875X0 1259 (X = Co, 


Ni). 


Composé Position du pic (eV) 


7.28 


11.45 


12.40 


Au-delà de 13 eV, l’amplitude de a(w) continue à diminuer et stabilise à la fin du spectre avec 
des valeurs égales à 1.97 (10*. cm") pour MgsS et 4.6 (10°. cm—') pour les deux composés 


Mgos75CO001258 et Mgo 875Ni0.1259- 
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Figure V. 9. (a) Conductivité optique o(o) et (b) coefficient d'absorption a(w) des composés Mes et 
Mgos875X0.1259 (X = Co, Ni), calculés par l'approximation TB-mBJ. 
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Conclusion générale 
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Conclusion générale : 


Dans cette thèse, nous avons étudié trois types de composés : 


1- Les composés Mg,,Cr,S (x= 0.125, 0.25, 0.50, 0.75) dans la structure rock salt 
ferromagnétique. Les propriétés calculées sont de types structurales, électroniques, élastiques, 


optiques et magnétiques. 


2- Les composés Mgs575X0125S (X= Cd, Zn) dans la structure rock salt nom magnétique. 
Les propriétés calculées sont de types structurales, électroniques, élastiques, optiques et 


thermoélectriques. 


3- Les composés Mg, 575X9125S (X —Co, Ni) dans la phase rock salt antiferromagnétique. 


Les propriétés calculées sont de types structuraux, électroniques et optiques. 


Les propriétés structurales des composés sont calculées par l’approximation WC-GGA, les 
propriétés électroniques, magnétiques et optiques par l’approximation TB-mBJ et les 


propriétés thermoélectriques par le code Boltztrap . 
Les résultats rencontrés ont montré : 


- que la valeur du paramètre de maille du composé binaire MgS pur concorde bien avec la 


valeur expérimentale publiée dans la littérature, 


- que les composés Mg,,Cr,S (x= 0.125, 0.25, 0.50, 0.75) sont des demi-métaux 
ferromagnétiques de polarisation en spin égale à 100% au niveau de Fermi et de moment 
magnétique total égal à 4 u4. Leurs propriétés élastiques ont montré qu’ils sont élastiquement 


stables et anisotropes. Les composés dopés aux concentrations x = 0,125, 0,25 et 0,50 sont 
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fragiles et le composé fortement dopé (x = 0,75) est ductile. Les propriétés optiques des 
composés ternaires ont mis en évidence une forte dispersion de la lumière dans la région 
ultraviolette, une nette différence a été observée entre les spectres optiques du composé 


binaire et des composés ternaires dans la région basse énergie. 


- que les composés Mg, 375X0125S (X= Cd, Zn) sont des semi-conducteurs non magnétiques 
caractérisés par une bande interdite large et sont mécaniquement stables. L'étude de leurs 
propriétés optiques a mis en évidence une absorption significative dans la région ultraviolette 
et l’étude de leurs propriétés thermoélectriques a montré qu’ils sont intéressants de ce point de 


vue à cause des valeurs importantes de leurs figures de mérite ZT. 


- que les composés Mgs375X91258 (X —Co, Ni) sont des semi-conducteurs à faible gap dans le 
canal des spins majoritaires et dans le canal des spins minoritaires et sont stables dans la 
phase antiferromagnétique, leurs propriétés optiques indiquent une réponse importante dans la 


région ultraviolette du spectre optique et une faible réponse dans la région visible. 


Les résultats obtenus de l’étude des composés réalisés par le dopage du Mg par les éléments 
magnétiques Cr, Cd, Zn et Co, indiquent que ces derniers présentent de bonnes 
caractéristiques qui leur permettent d’etre utilisés dans les dispositifs spintroniques, 


optoélectroniques et photovoltaïques. 
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